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FOREWORD 


This  document  Is  the  final  report  on  Contract  AF  19  (628)  -  2799,  It 
describes  an  experimental  program  devoted  to  the  study  of  the  particle  popu¬ 
lation  trapped  in  the  Earth's  magnetic  field  and  injected  into  it  by  either 
natural  cause  or  by  high  altitude  nuclear  bursts.  Instruments  were  constructed 
in  this  program  according  to  the  designs  developed  under  two  other  Air  Force 
contracts.  The  proton-alpha  detector  and  the  proton  spectrometer  were 
developed  under  Contract  AF  19  (604)  -  7347,  and  the  electron  and  gamma 
detector  was  developed  under  Contract  AF  19  (628)  -  252.  The  description 
of  these  instruments  is  abstracted  from  the  final  report  on  these  programs. 

Instruments  identical  to  the  ones  constructed  were  flown  on  an  Air 
Force  Satellite,  1962  ^k,  on  26  October  1962,  under  Contracts  AF  19  (628)  - 
252  and  AF  19  (604)  -  7347. 

The  analysis  of  the  data  obtained  by  this  instrumentation  has  been 
carried  out  as  part  of  this  program,  and  the  procedures  are  here  given. 

The  results  of  this  analysis  have  been  presented  to  the  scientific  community 
in  several  publications,  of  which  the  most  recent  are  enclosed. 


ABSTRACT 


This  document  describes  an  experimental  research  program  to  measure 
the  intensity,  angular  distribution,  and  energy  spectra  of  particles  trapped 
in  the  Earth's  magnetic  field.  Particular  attention  has  been  focused  on  the 
measurement  of  particles  artificially  injected  by  high  altitude  nuclear 
detonations  during  the  series  of  nuclear  tests  conducted  by  the  USSR  and 
the  USA  during  the  summer  and  fall  of  1962.  A  detailed  description  of  the 
instrumentation  and  calibration  procedures  are  given.  Measurements  obtained 
by  instrumentation  identical  to  the  one  here  described,  flown  on  an  Air 
Force  satellite  (1962  are  given.  The  method  of  data  handling  is  discussed 
and  machine  programs  for  data  reduction  are  described.  Results  from  these 
measurements  have  appeared  in  several  documents,  a  listing  of  which  is 
included  in  the  text.  A  summary  of  the  results  and  conclusions  is  given. 
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1.0  INTRODUCTION 


The  radiation  trapped  in  the  Earth's  magnetic  field  has  been  the  sub¬ 
ject  of  Intense  study  for  the  past  several  years.  Experiments  have  been 
prompted  by  the  intense  scientific  interest  in  this  important  geophysical 
phe’  menon  ais  well  as  by  the  practical  importance  of  its  effects.  Possible 
danger  to  manned  a.nd  unmanned  spacecraft  exposed  tc  the  trapped  radiation 
flux  has  been  a  primary  pxactical  coniideration.  Instruments  to  detect  and 
study  the  trapped  particle  environment  had  been  designed  and  flown  by 
several  groups  by  the  time  of  the  resumption  of  high  altitude  nuclear  testing 
in  19C2.  From  the  understanding  which  had  been  reached  on  the  behavior 
of  the  natural  radiation  belts,  it  was  rc'^Uzed  that  large  yield  high  altitude 
nuclear  detonations  would  significantly  alter  ''article  distribution  in  the 
belts  by  injecting  copious  numbers  of  the  charged  partiu*  produced  by 
the  detonation.  The  artificial  radiation  belts  thus  created  were  ...  ''act 
intense  enough  to  cause  damage  to  solar  cells  of  orbiting  spacecraft  and  ic 
present  a  real  hazard  for  manned  space  missions. 

An  important  difference  exists  between  natural  and  artificial  belts. 
Artificial  injections  result  in  a  sudden  increase  of  the  particle  intensities 
in  localized  regions  of  the  magnetosphere.  This  increase  can  be  so  large 
that  the  flux  of  artificially  produced  particles  is  orders  of  magnitude  greater 
than  the  natural  flux.  The  pitch  angle  distribution  and  energy  spectrum  of 
the  radiation  may  also  be  radically  changed.  With  time,  conditions  return 
to  the  natural  levels.  Depending  on  the  characteristics  of  the  injection  and 
on  the  processes  which  determine  the  removal  of  trapped  radiation,  the 
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characteristic  times  for  return  to  the  “nonnar  can  vary  by  large  factors. 
Estimates  of  these  times  prior  to  the  tests  varied  from  months  to  decades. 

The  purpose  of  the  research  program  described  in  this  rep^t  was  tc  stuoy 
the  nature  of  the  radiation,  its  angular  distribution,  its  energy  spectrum,  and 
its  temporal  and  spatial  variation.. 

Several  mechanisms  play  an  important  role  in  the  injection  of  trapped 
particles  produced  by  nuclear  detonation.  One  can  easily  understand  how  the 
decay  of  neutrons  or  of  radioactive  debris  can  result  in  trapping.  The 
relative  importance  of  the  various  processes  is,  however,  harder  to 
determine  and  depends  critically  on  the  depth  in  the  atmosphere  at  which 
the  detonation  occurs. 

The  low  atmospheric  density  at  the  altitude  at  which  Shot  Star  Fish 
Prime  was  detonated  permitted  the  debris  to  expand  to  great  distances. 

In  the  upward  direction  the  major  force  acting  to  restrict  this  expansion 
was  due  to  the  Interaction  of  the  charged  debris  with  the  geomagnetic 
field.  As  a  consequence,  this  expansion  continued  until  one  or  a  combina¬ 
tion  of  the  following  conditions  occurred.  The  debris  was  brought  to  rest, 
its  energy  having  been  converted  to  magnetic  energy  and  stored  in  the  com¬ 
pressed  field;  or  the  debris  was  guided  along  the  field  lines  of  the  mag¬ 
netic  bottUi  to  the  Southern  Hemisphere;  or  the  debris  de-lonized  and  the 
subsequent  trajectory  of  each  particle  was  determined  by  the  gravitational 
force  until  a  ^“decay  event  occurred;  or  the  expansion  became  t".rbuler.t 
with  the  debris  breaking  up  into  separated  blobs.  The  downward  expansion 
was  controlled  principally  by  the  exponentially  increasing  atmosphere. 

Several  ways  exist  to  permit  the  debris  to  expand  to  large  distances 
from  the  source.  Each  debris  fl-decay  event  results  in  the  Injection  of 
an  electron  into  the  magnetosphere.  What  happens  to  the  electron  depends 
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on  Its  energy,  the  oitltvxle  ot  injection,  and  thf^  pitch  angle  of  injection. 

For  each  altitude  electrons  injected  over  a  certain  portion  of  the  solid  anyle 
can  be  trapped.  The  contribution  by  this  p'ocesa  of  debris  decay  the 
artificial  belts  is  believed  to  have  been  predominant  ove”  other  processes 
and  in  particular  over  neutron  decay  injection  in  the  Star  Fish  event^^^ 

In  low  altitude  detonations  occurring,  for  instance,  at  30  km,  the  contri¬ 
bution  to  injection  by  the  netdron  decay  process  may  become  comparable 
or  predominate  over  the  debris  decay  Jnjectior  process.  Accurate  measure¬ 
ments  of  the  spectrum  and  angular  distribution  of  artificially  trapped  radia  * 
tion  could  permit  discrimlnaaon  between  the  different  injection  models. 

Of  great  importance  for  the  understanding  of  the  dynamics  of  the 
trapped  radiation  and  for  its  practical  consequences  is  the  determination 
of  the  removal  times  associated  with  the  artificial  radiation  belts.  The 
measurements  obtained  up  to  date  seem  to  indicate  a  rapid  initial  decay 
followed  by  s  slower  exponential  decay  with  characteristic  decay  times 
of  the  order  of  months.  It  would  also  appear  that  scattering  of  the 
trapped  particles  on  atmospheric  constituents  is  the  major  cause  for 
removal.^"  ^ 

The  results  of  the  program  here  desertbed  have  been  published  in 
various  forms.  In  this  report  a  brief  summary  of  results  and  of  the  major 
conclusions  are  given.  The  instrumentation  and  calibration  procedures 
are  described  in  detail. 
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2.0  INSTRUMENTATION 


The  instrumentation  constructed  under  this  contract  had  been  designed 
and  developed  under  other  contracts.  It  included  three  instruments: 

1.  A  proton  spectrometer  for  the  energy  region  of  6  to  70  Mev, 
developed  under  Contract  AF  19  (604)  -  7347.  The  detector 
and  pertinent  calibrations  are  described  in  Appendix  A. 

2.  A  proton  and  alpha  particle  detector  for  the  region  from  1  to 
10  Mev  also  developed  under  Contract  AF  19  (604)  -  7347  and 
described  in  Appendix  B. 

3.  An  electron  and  y  -ray  measuring  instrument  package  developed 
under  Contract  AF  19  (628)  -  252,  which  is  described  in  Appendix 
C.  This  instrument  contains  two  detectors  sensitive  to  electrons 
and  one  sensitive  to  y  -rays,  0.  05  to  5  Mev. 

The  three  instruments  were  constructed  under  this  program  to  replace 
identical  instruments  built  under  the  AF  19  (628)  -  252  and  AF  19  (604)  -  7347 
programs  and  flown  on  an  emergency  basis  on  an  Air  Force  satellite  given 
the  international  designation  1962  ^k,  flown  in  the  fail  of  1962. 

Part  of  this  contract  consists  of  the  analysis  of  the  data  obtained  by 
the  instrumentation  flown  on  this  satellite.  We  therefore  included  in  this 
report  a  description  of  the  principle  of  operation  and  characteristics  of  the 
three  instruments  mentioned  above. 
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3.  0  DESCRIPTION  OF  SATE  LUTE  MISSION 


In  the  fall  of  1962  an  Air  Force  satellite,  given  the  international 
designation  1962  ^k,  was  launched  into  orbit.  The  satellite  carried 
instrumentation  to  study  artificial  radiation  belts  created  by  the  injection 
of  electrons  from  nuclear  detonations  into  the  magnetosphere.  Under  the 
scientific  direction  of  Air  Force  Cambridge  Research  Laboratories,  a  num¬ 
ber  of  instruments  were  chosen  for  Inclusion  in  the  payload  The  instru¬ 
ments  described  in  the  preceding  chapter  were  part  of  the  payload. 

The  pertinent  parameters  for  the  satellite  orbit  at  early  times  were 
the  following:  apogee,  3000  n.m.;  perigee,  115  n.  m. ;  orbit  inclination, 
71°;  orbital  period,  147  minutes;  apsidal  rate,  0*8°  per  day.  Figure  1 
illustrates  these  orbital  parameters.  The  satellite  pitch  period  was 
127  seconds,  and  the  roll  period  was  53  seconds.  The  satellite  parameters 
were  well  suited  to  study  the  phenomena  of  interest.  For  example,  the 
high  eccentricity  of  the  orbit  enabled  the  acquisition  of  data  at  large  L 
values  near  the  equatorial  plane.  Slow  tumbling  rates,  along  with  the 
chosen  data  sampling  rate  of  about  one  per  second,  permitted  the  use  of 
instruments  which  could  measure  accurately  the  angular  distributions  of 
trapped  particles. 

"igure  2  shows  the  projection  of  an  early  orbit  on  a  dipole  field.  The 

(2) 

coordinates  are  the  usual  "natural"  B-L  coordinates  of  Mcllwain  .  The 
squares  and  ci.^cles  show  regions  of  data  acquisition  in  the  1. 2  and  1.  9 
L-shells,  shells  of  interest  with  respect  to  artificial  radiation  belts. 
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SATELLITE  ORBIT 


Figure  1 


LOCUS  OF  AN  ORBIT  IN  B-L  SPACE 


Th<i  instruments  comprising  the  total  payload  were  selected  to  yield 
redundant  measurements  in  the  energy  range  corresponding  to  fission 
electrcms.  and  to  include  narrow  acceptance  angle  detectors  as  well  as 
omni-directional  detectors,  instruments  obtained  their  power  from  a 
large  capacity  battery  pack  to  avoid  degradation  of  the  experiment  due  to 
possible  solar-cell  damage.  Data  could  be  acquired  both  by  direct  "real- 
time"  read-out  and  by  command  read-out  of  a  tape  recorder  capable  of 
storing  data  from  all  Instruments  for  a  time  longer  than  the  orbital  period. 

The  five  scintillation  detectors  comprising  the  AFCRL-ASE  portion 
of  the  payload  were  designated  Beta-1,  Beta-2,  Gamma,  P,  and  P-alpha. 
There  is  evidence  that  the  P-alpha  detector  malfunctioned  shortly  after 
launch.  The  most  salient  features  of  the  remaining  four  instruments  are 
summarized  in  Table  I. 

It  will  be  observed  that  detectors  Beta-1,  Beta- 2,  and  Gamma  have 
extremely  wide  fields  of  view  for  charged  particles;  (detector  Gamma  has 
essentially  a  4  ir  field  for  gamma  radiation).  These  instruments  were 
utilized  in  determining  the  gross  aspects  of  the  artificial  belts.  Detec¬ 
tor  P  has  an  extremely  narrov;  field  of  view,  and  is  essentially  a 
"directional  differential  intensity"  detector.  This  detector  was  used  in 
obtaining  pitch  angle  distributions.  Output  channels  most  used  for  the 
data  presented  later  in  the  report  are  marked  with  an  asterisk. 

Discrimination  among  types  of  particles  (electrons,  protons,  or 
gammas)  is  achieved  by  a  comparison  of  the  observed  counting  rates  and 
the  rates  of  ene'gy  deposition  In  several  channels,  and  by  roll -modulation 
observations.  In  artificial  radiation  belts  produced  by  high  altltiKie  nuclear 
detonations,  the  detected  particles  are  predominantly  electrons. 
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4.0  DATA  REDUCTION  TECHNIQUES 


The  telemetered  data  from  the  satellite  was  digitised  and  merged  with 
ephemerls  Information  by  Lockheed.  This  ephemetls  Information  included 
B-L  coordinates  derived  from  Mcllwaln's  model.  Programs  for  the  IBM  7090 
computer  were  developed  and  used  whenev^er  possible  In  the  data  reduction. 
These  programs  are  discussed  In  detail  In  Appendix  D. 

Digital  Computer  Processing.  Outputs  of  the  processing  system  Included: 

1.  Straight  listing  or  summary  of  telemetered  information. 

2.  B  -  L  matrix  with  detector  counting  rates  as  elements, 

3.  Pitch  angle  distribution  printouts. 

4. 1  Time  Correction 

It  was  discovered  that  the  universal  times  assigned  to  voltage 
points  on  the  digitized  tapes  were  in  error.  An  initial  assumption  was  then 
made*  to  ine  effect  that  the  time  error  was  constant  throughout  a  satellite 
orbit.  Then  two  methods  were  used  to  determine  time  corrections: 

1.  From  a  knowledge  of  the  satellite  orbit,  times  of  entry 
and  exit  from  the  Earth's  shadow  could  be  computed. 

This  information  was  correlated  with  readings  of  the 
satellite's  solar  sensors  to  produce  a  time  correction. 

2.  By  studying  the  instrument  readings,  it  was  possible  to 
determine  when  the  satellite  was  passing  through  radia* 
tlon  belts.  A  radiation  belt  has  a  constant  associated 


i 

i 

i 

L  value.  Since  ephemerls  quantity  L  was  computed  [ 

from  universal  time,  which  was  In  error,  successive 
passages  through  the  same  radiation  belt  had  apparent 
different  L  values.  By  adjusting  universal  tim^  (and 
ephemerls)  in  such  a  manner  as  to  yield  consistent  L 
values  on  each  passage  through  the  same  radiation 
belt,  a  tfme  correction  could  be  computed. 

Both  methods  yielded  consistent  results  for  any  satellite  ohase. 

We  have  therefore  used  in  data  analysis  universal  times  rorrected  by  using 
Method  1. 

4.  2  Correcting  Magnetometer  Data 

Figure  3  shows,  as  a  function  of  time,  a  typical  portion  of  the 
raw  data  from  detector  P  and  the  corresponding  calculated  value  of  the 
cosine  of  the  pitch  angle  6  ,  If  we  accept  these  pitch  angle  curves  at 
face  value,  we  would  obtain  inconsistent  distributions  from  the  left  and 
.ight  portions  about  a  maximum  or  minimum  of  the  cos  6  curve.  We  demand 
consistency  about  maxima  and  minima  in  our  analysis,  attributing  the 
apparent  an'-^ular  discrepancies  to  magnetometer  error.  This  amounts  to 
assuming  axial  symmetry  about  a  B-field  line. 

There  is  little  doubt  that  the  magnetometer  data  are  in  enor.  The 
angular  discrepancy  is  always  greater  in  low  B-field  regions.  Furthermore, 

If  we  plot  the  difference  in  absolute  values  between  the  magnetic  field 
strength  given  in  the  ephemerls  (calculated  from  the  Jensen  and  Cain  model) 
and  the  magnetic  field  strength  derived  from  the  magnetometers ,  we  obtain 
an  oscillatory  function  with  a  time  behavior  characteristic  of  the  pitch- 
and  roll  behavior  of  the  satellite.  Such  an  apparent  magnetometer  verror 
may  be  due  to  the  residual  magnetic  field  of  the  satellite  or  to  magnetometer 
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Figure  3 


cailbrfttion  drifts.  The  angular  discrepancy  does  not  seem  to  be  due  to 
timing  errors  since  in  the  sani3  orbit  the  angular  shift  required  to  bring  about 
data  consistency  is  iust  as  apt  to  be  positive  as  negative.  All  attempts  to 
discover  the  source  of  the  angular  discrepancy,  and  correct  for  it,  did  not 
prove  successful.  We  have,  therefore,  in  our  treatment  of  the  data,  shifted 
the  aspect  function  by  the  amount  required  to  bring  about  data  consistency. 
The  magnitude  of  such  a  shift  seldom  exceeds  6°  and,  as  stated,  tends  to 
be  larger  in  the  low  B-field  regions  of  the  orbit. 

A.SE  Document  ASE-845  gives  further  details  on  the  computer  programs, 
on  techniques  used  to  make  time  corrections,  and  on  attempts  to  make 
magnetometer  corrections. 
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5.0  RESULTS  AND  RECOMMENDATIONS 


After  reduction,  the  data  were  analyzed  to  achieve  the  following  aims: 

1.  To  establish  the  time  and  lc*cation  of  the  artificial  injections 
occuring  during  the  useful  life  of  the  spacecraft  instrumentation, 

2.  To  determine  the  angular  distributions  of  the  trapped  particles, 
and 

3.  To  study  the  temporal  and  spatial  intensity  variations. 

The  results  obtained  have  been  described  in  the  several  dociiments 
listed  below  to  which  ASE  scientists  have  contributed  under  this  program: 

1.  R.  Giacconi  and  F,  Paolini,  "Preliminary  Results  Obtained  by 
1401  Instrumentation  cn  the  Artificial  Radiation  Belts",  ASE 
Document  ASE-34 1 ,  1  January  1963.  SECRET 

2.  L.  Katz,  D.  Smart,  R.  Giacconi,  F.  Paolini,  and  R.  Talbot,  Jr., 
"Preliminary  Results  on  the  Artificial  Radiation  Belts  Since 
October  27,  1962,  as  Measured  by  Starad" ,  in  Symposium 
Proceedings:  Beta  Beta  and  Delta  Gamma  Programs,  6-7  March 
1963,  Armour  Research  Foundation  of  the  Illinois  Institute  of 
Technology,  Chicago,  Illinois,  DASA  1371,  DASA  Data  Center 
Special  Report  13,  May  1963.  SECRET-RESTRICTED  DATA. 

3.  L.  Katz,  D.  Smart,  F.  Paolini,  R.  Giacconi,  R.  Talbot,  Jr., 
"Measurements  on  Trapped  Particles  Injected  by  Nuclear  Detona¬ 
tions",  in  Symposium  Proceedings:  Trapf?ed  Radiation,  15-16  April 
1963,  Goddard  Space  Flight  Center.  Greenbelt,  Maryland, 
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DASA  1394,  DASA  Data  Center  Special  Report  16.  August  1963. 
SECRET. 

4.  L.  Katz,  D.  Smart,  F.  R.  Paollnl,  R.  Giacconl,  and  R.  Talbot,  Jr,, 
“Recent  Changes  in  the  Artificial  Radiation  Belts",  Presentaticxi  at 
the  44th  Annual  Meeting  of  the  American  Geophysical  Union, 
Washinc^on,  D.  C. .  18  April  1963, 

5.  L  Katz.  D.  Smart.  F.  R.  Paolini,  R.  Giaccont,  and  R.  Talbot.  Jr. , 
"Measurements  on  Trapped  Particles  Injected  by  Nuclear  Detona¬ 
tions",  Space  Research  IV.  Proceedings  of  the  Fourth  International 
Space  Science  Symposium,  Warsaw,  June  3-12,  1963,  North- 
Holland  Publishing  Company,  Amsterdam. 

6.  G.  A.  Kuck,  "Time  Behavior  of  Charged  Particles  Injscted  by 
1962  High  Altitude  Russian  Nuclear  Tests",  Thesis,  Air  Force 
Institute  of  Technology,  August  1964. 

Document  5  of  this  list  is  included  as  Appendix  E  of  this  report. 

In  summary,  the  results  of  the  program  were: 

1.  Detection  of  artificial  belts  due  to  Russian  High  Altitude  Nuclear  detona¬ 
tions  on  October  28  and  November  1. 

2.  Measurement  of  angular  distributions  and  lifetimes  of  the  trapped 
radiation  artificially  injected  by  the  Russian  detonations  irom  early 
times  up  to  45  days  after  the  detonation. 
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The  major  difficulties  encountered  in  the  programs  were: 


1.  Difficulty  of  Interpretation  of  the  data  due  to  the  broad  response 
characteristics  of  the  detectors  with  respect  to  different  species 
and  energies  of  particles. 

2.  Insufficient  angular  resolution  of  the  electron  detectors. 

3.  Incomplete  coverage  by  the  satellite  of  die  times  and  regions 
of  interest. 

4.  Insufficient  timing  and  aspect  accuracy. 

Most  of  these  difficulties  were  due  to  the  fact  that  the  instrumentation 
was  not  specifically  designed  for  monitoring  artificial  radiation  belts.  The 
short  times  involved  from  conception  to  flight  cf  1962  /3)r.  (about  i  month)  also 
did  not  permit  modification  of  the  detectors  or  adequate  discussion  and  im¬ 
plementation  of  the  timing  and  aspect  requirements.  Advanced  planning, 
resulting  in  stand-by  payloads  specifically  designed  to  detect  artificially 
injected  particles  and  ready  to  be  flown  for  such  missions  is  obviously  the 
answer.  This  policy  is  presently  being  imple  nented  by  the  Defense  Atomic 
Support  Agency. 

A  second  difficulty  which  as  yet  has  received  Inadequate  attention  is 
the  poor  coverage  which  one  can  obtain  from  polar  orbiting  satellites.  A 
polar  orbiter  affords  inspection  of  shells  with  a  wide  range  of  L  values 
but  poor  coverage  of  each  shell.  An  equatorial  satellite  with  identical 
instrumentation  with  appropriately  selected  orbit  could  permit  full  time 
coverage  from  zero  times. 
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APPENDIX  A 

"Proton  Spectrometer'’ 


Excerpt  from  G.  Davidson,  R.  Giacconi,  H.  Gur&Jcy,  et  al. , 

A  Research  Program  to  Investigate  the  Experimental  Problems 
of  Particle  Measurements  in  Space.  Final  Report  on  Contract  AF 
19  (604)  -  7347,  ASE  Document  ASE-400,  27  May  1963.  pp.  133-146. 
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3.  4  PROTON  SPECTROMETER 
3.4.1  Introduction 


The  knowledge  of  the  radiation  environment  in  which 
satellite  components  have  to  survive  is  of  particular  interest  for  the  design 
of  long-lived  systems.  Due  to  the  steepness  of  proton  and  electron  spectra 
in  the  Van  Allen  belt,  it  is  evident  that  most  of  the  energy  content  is  carried 
by  low -energy  particle  fluxes.  In  particular,  relatively  low-energy  protons 
(energies  greater  than  10  Mev)  can  traverse  considerable  thicknesses  of 
shielding  -  for  instance,  1  millimeter  of  aluminum  -  and  lose  their  total 
residual  energy  in  the  next  few  grams  of  material.  Damage  to  solar  cells 
by  such  protons  is  therefore  particularly  extensive. 

Such  data  must  be  known  for  satellites  which  spend  a  large  portion  of 

their  orbit  in  the  heart  of  the  inner  Van  Allen  belt  where  large  fluxes  of 

protons  in  this  energy  range  are  expected.  The  measurements  available 

to  date  have  been  perfcxmed  in  the  energy  range  above  60  Mev.  These 

data  indicate  a  steep  power  spectrum  up  to  320  Mev.  If  one  extrapolates 

these  results  to  the  10  Mev  region,  it  can  be  shown  that  10  Mev  protons 
3 

carry  about  10  times  more  energy  than  60  Mev  protons.  However,  such  an 
extrapolation  i-  not  known  to  be  valid,  and  it  is  therefore  of  great  importanc 
to  determine  the  differential  spectrum  of  protons  in  this  energy  range  to 
establish  the  existence  or  absence  of  a  "knee  "  in  the  e.nergy  spectrum 

The  following  sections  describe  an  instrument  which  ha.*;  been 
designed,  consiructeo  and  flown  primarily  for  the  purpose  of  measuring  the 
energy  spectrum  cf  protons  in  the  energy  region  6  to  100  Mev. 

3.4.2  Theory  of  Design 

The  guiding  philosophy  in  the  design  of  this  instrument 


was  to  achieve  the  utmost  simplicity  consistent  with  the  desired  objective 
of  measuring  the  spectrum  of  protons  in  the  energy  range  between  6  and 
100  Mev  in  the  heart  of  the  Inner  Van  Allen  belt.  The  spectrometer  is  a 
directional  instrument  which  measures  proton  fluxes  within  a  cone  +  5° 
from  its  axis,  and  therefore  is  capable  of  further  yielding  information  on 
pitch  angle  distributions. 

The  detector  consists  of  a  Thallium  activated  Csl  crystal  whose 

thickness  in  the  direction  of  the  incident  beam  is  equivalent  to  the  range 

-2 

of  90  Mev  protons  (11.  5  gr  cm  ).  This  sclntiilatcx  is  surrounded  by  a 
stainless  steel  shield  equivalent  in  thickness  to  about  100  Mev  protons 
except  in  a  narrow  opening  which  defines  the  field  of  view.  In  the  path 
of  the  incident  particles  and  in  front  of  the  CsI(Tl)  crystal  an  aluminum 

-2 

shield  equivalent  in  thickness  to  the  range  of  a  6  Mev  proton  (70  mg  cm  ) 
is  placed.  Figure  4  is  a  detailed  drawing  of  the  head  design. 

Protons  of  energy  greater  tnan  6  Mev  and  smaller  than  90  Mev  enter 
the  crystal  and  are  stopped  in  it.  The  crystal  therefore  has  a  response  which 
is  proportional  to  the  total  energy  of  the  proton  in  the  aforementioned  energy 
range. 

Protons  of  energy  greater  than  90  Mev  traverse  the  crystal  and  are  not 
stopped  in  it.  The  energy  loss  of  these  high-energy  particles  are  not 
linearly  proportional  to  their  total  energy,  out  equal  the  energy  loss  of 
some  lower  energy  protons  which  do  stop  in  the  crystal.  However,  in  anv 
particular  interval  of  energy  between  6  and  100  Mev,  the  amount  of  contamina¬ 
tion  due  to  such  high-energy  protons  will  be  less  than  10  per  cent  due  to  the 
steepness  of  the  proton  energy  spectrum.  Figure  5  shows  the  relationship 
between  the  Incident  proton  energy  and  the  pulse  height  obtained  from  the 
photomultiplier  viewing  the  CsI(T!)  crystal.  The  decrease  in  pulse  height 
with  increasing  energy  above  90  Mev  is  clearly  evident. 
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PROTON  SPECTROMETER  ASSEMBLY 
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This  detector  is  sensitive  elso  to  electrons  of  energy  greater  than 
‘10  kev,  which  can  penetrate  the  dlumlnurn  shield.  However,  since  the 
flux  of  high-energy  electrons  is  small  compared  to  the  proton  flux  of  the 
same  energy  in  the  inner  belt,  it  is  felt  that  the  electron  background  will 
not  give  rise  to  serious  difficulties  in  the  analysis  of  the  data  for  this  region 
of  space.  By  performing  pulse  height  analysis  of  the  scintillator  output 

it  is  therefore  possible  to  determine  the  differential  spectrum  of  protons 
between  6  3nd  100  Mev  to  reasonable  accuracy.  (  Energy  resolution  is 

about  15%  r\VHM,  and  readout  accuracy  about  jf  20%). 

3.  4.  3  Description  of  the  Instrument 

Figure  6  is  a  photograph  of  the  instrument.  It 

occupies  a  volume  of  6"  x  6"  x  12",  and  weighs  about  12  pounds.  The 
instrument  has  been  tested  to  operate  satisfactorily  over  a  temperature 
range  of  20°F  to  165°F,  and  can  withstand  the  vibration,  shock,  and 
acceleration  environment  specified  for  Midas  vehicles. 

Figure  7  gives  a  block  diagram  of  the  instrument.  The  scintillator 
is  viewed  by  a  ruggedized  photomuiiioiier.  The  average  dc  current  of  the 
photomultiplier  anode  is  monitored  in  addi  -on  to  the  singl(2  pulses  and  is 
therefore  a  measure  of  the  total  rate  of  energy  defxjsitlon  in  the  cr/stal. 

The  dc  current  is  fed  to  a  dc  amplifier  which  presents  the  information  in 
analog  form  between  0  and  5  volts  to  a  telemetering  point. 

Pulses  are  fed  Uirough  two  amplifiers  with  different  gains.  In  the  first, 
pulses  correspondina  to  energy  losses  in  the  scintillator  between  1  and  8  Mev 
are  amplified  to  give  pulse  heights  between  0.5  and  4.0  volts.  In  the 
second,  pulses  corresponding  to  energy  losses  between  20  Mev  and  70  Mev 
are  amplified  to  give  pulse  heights  between  1.0  and  3.5  volts.  The  outputs 
of  these  two  amplifiers  are  discriminated  over  eight  ranges  of  pulse  height 
corresponding  to  energy  losses  greater  than  1,  2,  4,  8,  20,  30,  >0  and 
70  Mev.  Figure  5  illustrates  the  relationship  between  incident  proton 
energy,  energy  loss,  anri  approf’irlate  discrimination  levels. 
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The  output  of  each  energy  channel  is  fed  to  a  logarithmic  rate  meter 
whose  voltage  is  proportional  to  the  logarithm  of  the  counting  rate  over  a 

dynamic  range  of  3  (xders  of  magnitude.  In  the  low-energy  channels,  this 

4 

range  extends  from  10  to  10  counts  per  second  while  in  high-energy 

3 

channels  it  extends  from  0.  1  to  10  counts  per  second.  The  output  of  the 
logarithmic  rate  meters  is  in  analog  form  between  0  and  5  volts.  The  different 
ranges  in  sensitivity  were  chosen  to  cover  the  expected  possible  variations 
in  the  measured  intensities  which  might  be  encountered  in  the  heart  of  the 
inner  Van  Allen  belt. 

Twelve  telemetering  points  are  needed  to  transmit  the  analog  infcx’mation. 
Additional  telemetering  points  are  used  to  monitor  the  photomultiplier  high 
voltage  and  the  transistcx  power  supply  voltage. 

C’rcuit  diagrams  i(x  the  various  elements  comprising  the  instrument 
are  identical  to  corresponding  elements  of  the  P-a  detector,  and  will  not 
be  duplicated  here.  The  power  consumption  of  the  total  instrument  is  about 
2  watts. 

3.4.4  Calibration  Procedures 

Because  of  time  limitations,  it  was  impossible  to  calibrate 
the  two  instruments  of  this  type  with  proton  accelerators. 

Two  acceleratOTs  would  have  been  necessary  -  a  low  energy  machine  (such 
as  the  8  Mev  ONR  Van  de  Graaff  at  M.  I.  T. )  and  a  high  energy  machine 
(such  as  the  200  Mev  Harvard  Cyclotron)  -  for  calibration  of  the  low  and 
high  energy  channels  respectively. 

_  „  210 

ior  tnis  reason,  a  Po  alpha  particle  source  was  used  for  energy 
calibration.  The  source  was  in  the  shape  of  a  very  thin  Po  disk,  one-eighth 
inch  in  diameter,  covered  with  a  thin  protective  film  of  evaporated  gold, 
ail  mounted  on  an  aluminum  planchette.  The  energy  of  the  alpha  particles 
penetrating  the  gold  was  measured  to  be  5.  2  Mev  with  a  0.  2  Mev  half¬ 
maximum  width  by  the  suppliers  of  the  source,  (The  spread  is  caused  by 
energy  losses  in  the  source  and  the  protective  gold  foil.)  This  energy 
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spreac  is  well  within  the  resolution  of  scintillation  detectors,  and  the  source 
wcs  therefore  perfectly  adequate  for  our  purpose- 

calibration  of  the  eight  flux  (counting  rate)  channels  of  each  instrument 
was  achieved  as  follows.  The  ohotomultiplier  high  voltage  was  preset.  The 
gain  of  the  amplifier  feeding  the  four  low  energy  channels  was  then  adjusted 
until  the  amplifier  pulses  produced  by  the  alpha  source  in  contact  with  the 
scintillates  were  1  65  vclts  high.  Since  a  5.  2  Mev  alpha  particle  gives 
the  same  pulse  height  in  CsI(Tl)  as  a  3.3  Mev  proton,  this  set  an 
amplifier  scale  factor  of  0.  50  volts /Mev  foe  protons.  The  photomultiplier 
was  then  disconnected  from  the  amplifier  input,  and  a  pulse  generator  with  a 
network  to  vield  the  same  pulse  shape  as  was  obtained  from  the  photomultiplier 
was  substitutea.  The  discriminaters  of  the  four  low  energy  channels  were 
then  set  to  fire  for  amplifier  output  pulse  heights  of  0.  5,  1.0,  2,  0,  and  4.  0 
volts,  corresponding  therefore  to  protons  which  give  energy  losses  in  the 
scintillator  of  1,2,  4,  and  8  Mev  respectively.  The  relationship  between 
incident  proton  energy  and  energy  loss  in  the  scintillator  (indicated  in 
Figure  5  was  calculated  from  proton  range  -  energy  relationships  in  the 
usual  fashion,  using  the  Univ.  of  Calif.  Range-Energy  Tables,  UCRL-2301. 

The  gain  of  the  amplifier  feeding  the  four  high  energy  channels  was 
adjusted  until  the  gain,  relative  to  the  low  channel  amplifier,  was  just  0.  10. 
This  was  accomplished  in  practice  by  trimming  this  amplifier  until  a  0.  5  volt 
pulse  height  output  was  obtained  when  the  output  from  the  low  energy  channel 
amplifier  was  5.  0  volts;  the  test  pulser  was  used  to  furnish  the  input  pulses. 
This  procedure  thus  set  the  high  energy  channel  amplifier  scale  factor  to 
0.  05  volts/Mev  for  protons.  The  discriminators  of  the  four  high  energy 
channels  were  then  set  to  fire  for  amplifier  output  pulse  heights  of  1.0, 

1.5,  2.  5,  and  3.  5  volts,  corresponding  to  protons  which  give  energy  losses 
in  the  scintillator  of  20,  30,  50,  and  70  Mev  respectively.  Since  both 
aripiifiers  were  tested  to  be  linear  within  a  few  per  cent  over  output  pulse 
heights  from  0.  2  to  5.  0  volts,  this  procedure  resulted  in  a  fairly  accurate 
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calibration  of  the  high  energy  channels  despite  the  indirectness 
of  the  technique. 

Calibration  of  the  energy  deposition  channel  was  achieved  as  follows. 
The  instrument  was  placed  in  a  standard  location  in  front  of  an  opening  in  a 
lead  house  containing  a  10.  milllcurie  Co  source.  The  rate  of  energy 
deposition  in  a  standard  shape  CsI(Tl)  crystal  had  previously  been  determined 
for  this  location  In  the  following  way.  The  crystal  was  mounted  on  a  photo¬ 
multiplier,  run  at  a  predetermined  negative  high  voltage  and  the  charge 

22 

deposited  by  the  Na  0.  51  Mev  annhilation  photopeak  was  measured  in  the 

conventional  way.  This  led  to  a  "charge:  sensitivity  constant",  giving  the 

number  of  coulombs  deposited  per  ev  of  energy  deposition  by  electrons.  This 

number,  of  course,  is  also  identical  in  magnitude  to  a  "current  sensitivity 

constant"  giving  the  photomultiplier  current  (coulombs /second)  produced 

per  ev/second  rate  of  energy  deposition.  The  PM  was  then  rewired  for  anode 

60 

dc  current  measurement,  and  placed  at  the  standard  location  before  the  Co 
source.  The  PM  was  run  at  the  same  high  voltage  as  before,  and  the  anode 
current  was  measured.  Since  the  "current  sensitivity  factor"  was  known, 
this  led  directly  to  the  rate  of  energy  deposition  in  the  standard  CsI(Tl) 
cr>'stal  at  the  standard  location  before  the  source.  Geometrical  and  shielding 
corrections  were  then  calculated  for  the  CsI(Tl)  scintillator  and  shielding 
proper  to  the  proton  spectrometer  to  be  calibrated.  These  were  applied  to 
derive  a  new  source  strength,  in  ev/sec,  peculiar  to  this  detector  placed  in 
the  standard  location.  The  signal  output  from  this  instrument's  current 
meter  was  then  recorded.  This  procedure  calibrated  one  point  for  this  channel, 
in  an  absolute  fashion.  To  achieve  calibration  over  the  entire  dynamic 
range  of  this  channel,  use  was  made  of  the  fact  that  photomultiplier  anode 
current  (within  space  charge  limitations)  is  linearly  proportional  to  the  rate 
of  energy  deposition  in  the  scintillator.  An  intermediate  calibration  curve, 
of  current  meter  signal  output  versus  input  current  (generated  by  a  calibrated 
laboratory  current  source)  was  first  plotted.  The  recorded  photomultiplier 
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signal  obtained  from  the  Co^^  source  led  directly  to  the  needed  proportionality 
constant  between  PM  anode  current  and  rate  of  energy  deposition.  The  current 
scale  of  the  aforementioned  intermediate  calibration  curve  was  then  relabbd  in 
terms  of  the  rate  of  energy  deposition,  resulting  in  an  absolute  calibration 
curve.  The  curve  for  a  typical  instrument  is  shown  in  Figure  8  . 

The  logarithmic  count  rate  meters  were  calibrated  by  the  use  of  the 
pulser  utilized  in  the  energy  scale  (discriminator)  calibrations.  Output  signal 
versus  the  pulse  repetition  rate  input  was  plotted  for  each  of  the  eight  flux 
channels  of  the  instrument.  Figure  9  shows  a  typical  calibration  curve. 

3.4.5  Summary  of  Experimental  Results 

Up  to  the  time  of  this  report,  one  protcn  spectrometer  has 
flown,  and  that  one  on  a  satellite  given  the  international  designation  1962/3it. 

The  instrument  operated  satisfactorily.  Analysis  of  the  data  obtained  has  not 
yet  been  completed.  Preliminary  results  were  reported  at  the  NASA-DASA  Sym¬ 
posium  on  Trapped  Radiation,  held  at  Goddard  Space  Flight  Center  in  Greenbelt, 
Maryland,  on  April  15  and  16;  and  again  at  the  44th  Annual  Meeting  of  the 
American  Geophysical  Union  held  in  Washington,  D.C.  ,  cn  April  17  through  20. 

Probably  the  most  interesting  data  obtained  from  this  instrument  resulted 
from  its  capabilities  as  an  electron  spectrometer,  rather  than  as  a  proton  spec¬ 
trometer.  In  this  application,  since  the  electron  to  proi^n  pulse  height  ratio 

2 

is  very  close  to  unity  in  CsI(Tl),  and  since  the  70  mg/cm  aluminum  shield  in 
front  of  the  scintillator  results  in  negligible  energy  loss  for  high  „nergy  electrons, 
the  Urst  five  disv-riminator  levels  c  orrespond  verv'  nearly  to  electrons  of  incident 
energy  1,  2,  1,  H,  and  20  Mev.  (The  last  three  discriminators  will  not  resfxmd 
to  electrons  since  the  maximum  energy  loss  an  electron  can  experience  in  the 
CsI(Tl)  is  abou-  25  Mev.)  Durlno  flight,  significant  electron  counting  rates  were 
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observed  in  the  first  three  channels.  Highest  fluxes  occurred  when  the  satellite 
passed  through  the  L  -  1. 83  shell,  containing  an  artificial  radiation  belt.  This 
L  region  is  also  the  slot  between  the  natural  inner  and  outer  radiation  belts, 
known  to  be  low  in  proton  population.  Correlation  with  outputs  of  other  instru¬ 
ments  flown  established  conclusively  the  identification  of  the  particles  beyond 
L  =  1.  8  as  electrons.  Protons  were  detected  by  the  instrument  in  the  region 
of  the  inner  belt  as  expected,  but  data  analysis  for  this  region  of  the  magneto¬ 
sphere  is  incomplete. 
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Excerpt  from  G.  Davidson,  R.  Giacconi,  H.  Gursky,  et  al. , 
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of  Particle  Measurements  In  Space,  Final  Report  On 
Contract  AF  19  (604)  -  7347,  ASE  Document  ASE-400, 

27  May  1963,  pp.  58  -  89. 


3.  1.  2  Instrumentation 


The  use  of  an  artificial  satellite  as  a  vehicle  imposed 
particular  requirements  on  the  weight,  reliability,  and  ruggedness  of  the 
instrumentation  necessary  to  perform  the  experiment.  Classical  solutions 
for  mass  aiscrimination  such  as  magnetic  separation  had  to  be  discarded 
owing  to  weight  considerrtion.  In  order  to  exploit  the  benefits  of  long 
available  times  of  observation,  the  power  consumption  requirements  were 
kept  to  a  minimum.  Particular  problems  existed  in  the  data  transmission. 

The  experimental  arrangements  which  are  discussed  in  the  following 
sections  were  designed  to  satisfy  these  requirements. 

Emphasis  was  placed  on  simplicity  of  construction,  reliability  of 
operation,  light  weight,  and  low  power  requirements.  A  particular  effort 
was  spent  to  reduce  to  a  minimum  the  requirements  on  the  data-transmltting 
equipment. 

The  instrumentation  required  for  the  proposed  experiment  was  designed 
to  perform  any  or  all  of  the  following  functions: 

A.  Detect  protons,  alpha  particles,  and  other  possible 
heavT  components  in  the  Mev  energy  range, 

B.  Distinguish  between  these  particles. 

C.  Measure  their  energy. 

The  three  functions  are  not  independent  from  each  other.  In  fact,  to 
discriminate  betv/een  different  particles,  a  measurement  of  ..he  energy 
together  with  another  independent  quantity,  such  as  range,  or  rate  of 
energy  loss  is  required.  Therefore,  rather  than  discussing  single 
components  which  could  be  used  to  perform  one  particular  function,  we 
prefer  to  discuss  complete  systems  which  can  accomplish  all  three  functions. 
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Each  such  system  is  built  around  a  basic  component,  the  detector. 

The  characteristics  of  the  detector  and  the  type  of  information  it 

can  furnish  dictated  the  choice  for  additional  components  which  had 

to  be  used  to  perform  the  experiment.  For  the  final  design  of  the  experiment, 

various  possible  systems  built  around  adaptions  of  available  detectors 

were  studied  in  detail  and  a  scintillator  detector  was  selected  as  the  basic 

detector. 

3.  1.  3  Principle  of  Operation 

Inorganic  scintillators  such  as  thallium-activated  Csl 
will  respond  almost  linearly  to  the  energy  loss  within  their  volume.  By 
using  more  than  one  scintillation  detector,  one  could,  therefore,  measure 
both  the  specific  energy  loss  and  the  total  energy  of  the  particle.  In 
Figure  10a,  a  schematic  arrangement  is  shown.  The  anti-coincidence 
detector,  S  ,  assures  that  the  particles  stop  in  and,  therefore,  that 
the  pulse  height  in  is  a  measure  of  the  total  energy,  E.  The  pulse 
height  in  is  a  measure  of  the  specific  energy  loss  AE/Ax.  From 
these  two  measurements,  the  type  and  energy  of  the  particle  can  be 
determined.  This  classical  arrangement,  however,  is  easily  adapted 
only  to  the  detection  of  particles  of  relatively  high  energy.  For  low- 
energy  particles,  the  thickness  of  S,  must  be  made  smaller  than  the 
residual  range  of  the  particle.  In  the  case  of  protons  or  alpha  particles 
of  energies  of  a  few  Mev,  the  residual  ranges  are  of  the  order  of  microns, 
as  shown  in  Figure  11. 

The  practical  difficulties  involved  with  a  thin  detector  led  to  a 
consideration  of  other  possible  schemes.  A  different  approach  which  was 
tire  one  finally  adopted  consists  in  the  measurement  of  maximum  range 
(rali.er  than  AE/dx)  and  energy,  ^nis  arrangement  is  the  one  shown 
schematically  in  Figure  10b.  The  counter  which  previously  was  used 
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to  measure  AE/Ax  is  eliminated.  One  can  only  measure  the  total 
energy^  £,  lost  in  S^.  Still,  it  can  be  shown  that  such  a  detector 
would,  in  fact,  be  able  for  given  energy  regions,  to  distinguish 
between  different  kinds  of  particles  besides  measuring  their  energy. 

Let  us  assume,  for  instance,  that  the  scintillator  $2  has  a  given 
thickness,  R.  The  maximum  range  that  a  particle  stopping  in  $2 
can  have  is  (neglecting  straggling)  also  R,  Different  particles 
having  the  same  range  have  different  energies,  as  evident  by  inspection 
of  Figure  11.  For  Instance,  the  maximum  energy  of  a  proton  stopping 
in  .  2  millimeters  of  CsI(Tl)  is  about  5  Mev,  while  the  maximum  energy 
of  an  alpha  particle  stopping  in  this  thickness  is  20  Mev.  Therefore, 
no  pulse  from  corresponding  to  an  energy  loss  of  mere  than  5  Mev  can 
be  given  by  a  proton.  Consequently,  with  this  particular  thickness  of 
scintillator,  we  are  able  to  separate  alpha  particles  of  energies  between 
5  and  20  Mev  from  protons  of  any  energy,  and  are  also  able  to  measure 
the  a-particle  energy  spectrum  *n  tliis  region.  With  other  similar  detectors 
having  different  scintillator  thicknesses,  other  regions  of  the  spectrum 
could  be  selected.  It  should  be  noted  that  the  higher  Z  particles  could 
not  be  separated  from  the  alpha  component  in  this  fashion,  so  that  a  small 
contamination  of  high  Z  would  not  be  resolved.  The  same  method  could, 
however,  be  used  to  separate  high  Z  components  from  alpha  particles  in 
given  spectral  regions.  The  design  of  this  particular  type  of  system  was 
adopted  for  the  flight  instrumentation. 

3.1.4  Rejection  of  Spurious  Events 

The  available  measoreroents  on  intensity  and  composition 

of  radiation  in  the  V'an  Allen  belts  of  radiation  indicate  omni-directional  in- 

9  2 

tensities  in  the  core  of  the  inner  layer  of  10  /cm  see.  for  electrons  with 

4  2 

energies  greater  than  20  kev  and  2x10  /cm  sec.  for  protons  with  energy 
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greater  than  about  40  Mev.  While  the  electronics  associated  with  the 
detector  can  ver>  effectively  distinguish  between  the  pulse  heights 
produced  by  the  electrons  and  heavier  particles,  it  was  felt  that  the 
number  of  electrons  present  is  so  high  that  they  might  produce  a  dc 
background,  which  would,  in  effect,  jam  the  counter.  One  could  improve 
this  situation  by  reducing  the  solid  angle  viewed  by  the  detector,  but  in 
this  way,  the  fluxes  of  ali  types  of  particles  would  be  reduced  by  the  same 
amount.  We  used  a  different  approach  to  this  problem  consisting  of  the 
use  of  a  magnetic  separator  capable  of  reducing  the  electron  fluxes  by 
orders  of  magnitude  without  appreciably  reducing  the  proton  and  alpha 

fluxes.  A  magnetic  field  of  about  2000  gauss  ever  a  region  of  6. 

2 

cm  was  used  to  deflect  particles  with  momenta  up  to  10  Mev/c 
by  30°.  Electrons  up  to  10  Mev  energies  could  thereby  be  re¬ 
moved  from  the  accepted  beam,  while  protons  and  alphas  of  the  same 
energy  would  be  deflected  by  angles  of  the  order  of  1°.  The  advantage 
of  this  method  is  due  to  the  fact  that  proton  and  alpha  fluxes  are  not 
attenuated,  while  the  electron  flux  is.  The  magnetic  field  required  was 
obtained  by  using  small  permanent  magnets.  A  diagram  of  the  detector 
assembly  is  shov/n  in  Figure  12. 

Care  was  exercised  in  the  design  of  the  experiment  to  discriminate 
against  other  possible  sources  of  background  as,  for  instance,  nuclear 
interactions  of  low  and  high-energy  particles  in  the  detector.  3y  using 
a  thin  detector  we  reduced  the  probability  of  occurrence  of  such  interactior^ 
to  relatively  low  values. 

3.  1.  5  Flight  Instrumentation 

Figure  13  Is  a  photograph  of  the  flight  instrumentallori 
which  was  constructed  to  Implement  this  experiment. 

A  listing  of  the  design  specifications  is  given  below. 
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ALPHA- PROTON  PARTICLE  DETECTOR 


Power  Requirements; 
Weight: 

Detector: 

Field  Strength: 
Shielding: 

Geometric  Factor: 
Max.  Counting  Rate. 
Counting  Life: 
Outputs: 


Environmental 
Operating  Limits: 


1  watt,  28  +VDC. 

10.  8  lbs. 

6199  photomultiplier  with  Csl  and  plastic  scintillator 
sandwich. 


B  500  gauss. 

1.5  m  gr/cm^ 

-3  2 

Ail=  1.  5  X  10  cm  stcr. 

7  2 

10  protons/cm  sec  ster. 

Indefinite. 


There  are  eight  outputs  fed  to  th.e  cable  connector. 
These  consist  of  three  log  count  rate  meters  whose 
dc  output  is  a  linear  function  of  the  logarithm  of 
their  counting  rate.  The  discriminator  circuits  in 
these  channels  are  set  to  respond  to  pulses  of  proton 
and  alpha  particle  stopping  in  the  Csl.  The  energy 
loss  required  to  trigger  the  discriminators  are  as 
follows: 

Ch,  Pi  >  1  Mev 

Ch.  P2  2  Mev 

Ch.  P3  i  3.  5  Mev 


In  addition  to  this,  there  is  the  output  from  a  cyclic 
gate  (Ch.  1)  which  internally  commutates  the  outputs 
from  four  other  log  count  rate  meters  and  displays 
each  for  about  four  seconds.  These  channels  respond 
to  particles  having  the  following  energies: 


Ch.  o  la 
Cn.  a  2a 
Ch.  hp  2a 
Ch.  hp  la 


5.  5  -  16  Mev  alphas 
10  -  15  Mev  alphas 
10-50  Mev  protons 
15  -  Mev  protons 


In  conjunction  with  this  output,  and  having  four 
separate  outputs,  are  four  digital  readouts;  Ch  a  lb, 
Ch,  o  2b,  Ch,  hp  2b,  Ch.  hp  lb.  each  operated  by 
an  analog  staircase  and  covering  the  same  energy 
ranges  as  the  above  log  count  rate  meters 


Temperature: 

Acceleration; 

Shock;: 

Vibration; 


+  20°r  to  +165°F 
SOG's 

50  G‘s  in  7  X  10  "  sec. 
25  G’s.  30  -  3000  cps. 


A  block  diagram  of  the  electronics  is  shown  iii  Figure  14.  (Note: 

Some  of  the  energy  ranges  indicated  differ  from  those  of  the  previous  table; 
the  table  indicates  the  later  design  values). 

Figures  15,  16,  17,  18,  19,  20,  and  21,  show  the  wiring  of  the 
individual  Modules  A  Through  G. 

A  particle  in  the  field  of  view  of  the  detector  is  deflected  by  the 
static  magnetic  field.  Electrons  are  deflected  enough  so  that  approximately 
only  1%  of  the  electrons  at  1  Mev  can  reach  the  detector.  Proton  and  alpha 
particles  are  deflected  almost  inappreciably. 

A  particle  impinging  on  the  detector  can  penetrate  one  or  both  of  the 

scintillators.  A  m.ixed  light  pulse  from  the  Pilot  B  scintillator  with  about 

“8  ’•7 

5x10  sec  decay  time  and  fi  .*m  the  CsI(Tl)  with  about  3  x  10  sec  decay 

time  will  be  detected  and  amplified  by  the  photomultiplier.  Pulses  are 

monitored  at  the  last  dynode  and  at  the  anode.  The  dynode  pulse  is  clipped 

by  use  of  a  shorted  cable  which  eliminates  the  slow  component. 

The  anode  pulse  is  amplified  and  fed  to  a  linear  gate  (Module  C2)  which 
is  triggered  by  the  identical  pulse  delayed  by  .  2  x  10  ^  sec.  This  allows  only 
the  slow  portion  of  the  pulse  to  be  transmitted  to  an  antigate  (Module  B)  which 
IS  triggered  by  the  fast  pulse.  If  no  "fast  pulse"  is  present  the  slow  component 
from  the  CsI(Tl)  is  fed  to  a  5  channel  analyzer  (Modules  D1  through  D5).  If  a 
fast  pulse  is  present  no  output  is  present  from  Module  B. 

On  the  other  hand  both  pulses  are  also  fed  to  a  linear  gate  (Module  Cl). 

The  fast  pulse  is  delayed  by  .  5  x  10  ^  sec  to  allow  coincidence  between  it 
and  the  slow  portion  of  the  pulse  from  (Module  C2)  the  previous  linear  gate  which 
furnishes  the  slow  portion  and  constitutes  the  trigger.  If  coincidence  exists 
the  last  pulse  is  analyzed  on  2  channels  (Modules  D6  and  D7).  The  meaning 
of  these  operations  from  the  point  of  view  of  the  measurement  is  the  following. 
The  pulsi  from  a  slow  particle  stopping  in  the  first  scintillator  Csl(T^)  is 
analyzed  by  D1  to  D5.  The  pulse  height  corresponds  to  its  total  energy. 
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BLOCK  DIAGRAM  OF  ELECTRONICS 


MODULE  A 


(PIN  M) 


MODULE  "B"  DISCR.  FAST  UNIVIBRATOR  8  LINEAR  GATE 


1U21 


MODULE  E 


2N706  2N706 


The  pulse  from  a  fast  particle  which  enters  both  scintillators  is 
analyzed  by  D6  and  The  level  of  discrimination  in  the  linear  gate 
is  such  that  only  particles  which  exhibit  a  large  enough  energy  loss 
in  l!»e  CsI(Tl)  are  accepted.  This  corresponds  to  a  AE/Ax  requirement. 

The  level  is  such  that  only  particles  with  energies  less  than  50  Mev  are 
accepted.  These  particles  will  stop  In  the  Pilot  B  scintillator  if  they  are 
protons  of  less  than  32  Mev  or  electrons  of  less  than  2  Mev.  By  requiring 
.  5  Mev  energy  loss  in  the  first  crystal  and  greater  than  10  and  15  Mev 
energy  loss  in  the  second,  we  can  eliminate  electrons  and  measure  the 
energy  of  the  protons  in  two  bands  of  10  to  50  and  15  to  30  Mev.  The 
outputs  of  all  discriminators  are  fed  to  count  rate  meter  which  furnishes 
a  voltage  approximately  proportional  to  the  logarithm  of  fee  counting  rats 
(Module  E).  Atypical  callbratior  '•urve  is  shown  in  ?igure  22.  Also 
on  the  four  highest  energy'  channels  the  output  is  fed  to  octal  staircase 
scalers  (Module  F)  to  accurate  chanting  at  very  low  fluxes.  The 
output  of  these  scalers  were  multiplexed  within  the  experiment  by  means 
of  a  cyclic  gate  (Module  G). 

Additional  electronics  included  a  high  voltage  supply  of  the  Ccckcroft 
Walton  type  and  a  dc  to  dc  converter  and  voltage  regulator  shown  in 
Figure  23.  Figure  24  shows  the  wiring  diagram  of  a  test  panel  for 
testing  Modules  A  through  G. 

The  output  from  the  channels  were  adjusted  to  measure  particles  in 
the  follov.’ing  ranges  of  energy^  given  on  page  47, 
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28V  TO  6.75V  DC-DC  CONVERTER 


WIRING  OF  PANEL  FOR  TEST  OF  MODULES 


3 .  L  6  Calibrations 

The  prototype  flight  instrument  was  subjected  to 
extensive  testing.  Exposure  at  the  High  Voltage  Engineering  Van  de 
Graaff  electron  accelerator  was  used  to  study  the  response  of  the 
instrument  to  very  high  electron  and  gamma  ray  fluxes.  Effects  such 
as  pile-up  and  fatigue  were  studied  to  assure  that  conservative  criteria 
had  iDeen  used  in  the  design  of  the  experiment. 

Both  in-house  calibrations  with  alpha  particle  sources  and  calibra¬ 
tion  at  the  14  Mev  M.I.T.  Cyclotron  were  used  to  set  the  proper  energy 
discrimination  levels.  The  protons  and  alpha  particles  beam  was  brought 
outside  through  a  6.4  x  10  gr/cm  Au  foil.  The  energy  of  the  particles 
was  then  attenuated  by  use  of  A1  foils  .001  inches  thick  in  various 
amounts.  Protons  with  energy  loss  as  small  as  .  5  Mev  in  the  scintilla¬ 
tors  were  detectea.  A  point  calibration  was  performed  to  determine 
linearity  of  the  system  to  a  given  energy  lost  in  the  CsI(Tl).  Figure  25 
shows  a  tvyical  curve  obtained  in  such  a  test.  The  vertical  error  bars  are 
estimated  on  the  basis  of  inherent  statistical  accuracy  of  the  detector.  The 
horizontal  error  bars  are  due  to  uncertainty  inherent  in  the  method  used  to 
degrade  the  energy.  Effects  of  straggling  at  very  low  energies  are  predomi¬ 
nant.  The  efficiency  for  measuring  alpha  particles  was  determined  and  a 
resolution  curve  obtained  with  alpha  particles  showed  an  energy  resolution 
of  about  20  per  cent.  The  rejection  efficiency  of  the  anti  coincidence  was 
measured  and  found  to  be  identical,  within  the  measurement  error,  to  IOC 
per  cent.  These  calibrations  were  repeated  for  each  flight  unit  to  assure 
reproducibility.  It  was  found  that  the  individual  differences  in  geometry  and 
crystal  efficiencies  v/ere  negligible. 

In  the  field,  calibration  was  carried  out  by  espsciaiiy  mounted  radioactive 
alpha  particle  sources  which  could  be  brought  into  close  proximity  of  the 
CsI(Tl)  crystal. 


61 


OUTPUT  PULSE  (VOLTS) 


ENERGY  CALIBRATION  CURVE 
AT  lyl.l.T.  CYCLOTRON 


ESTIMATED  ENERGY  LOST  IN  Csl  (Tt)(M«v) 

62 


3.1.7  Operational  Information 


Seven  complete  units  of  the  paroton-alpha  experiiner*ts 
(AFCRL  designation  CRM -28)  were  built.  The  first  unit  was  delivered  on 
1  February  1961. 

Unit  ASE-P-32-0  was  used  in  destructive  environmental  testing. 

Unit  ASE-P-32-1  was  flown  on  the  Atlas  Pod  08,  at  Cape  Canaveral 
on  19  December  1961. 

Unit  ASE-P-32-4  was  flown  on  an  Air  Force  Satellite  launched  on 
14  December  1962  at  VAFB  (PMR).  The  satellite  did  not  achieve  orbit. 

Unit  ASE-P-32-6  was  (under  the  direction  of  the  scientific  monitcx) 
modified  to  fly  on  the  1962  3k  satellite.  This  unit  was  working  properly 
up  to  a  few  hours  before  launch  when  last  calibrated.  The  unit  did  not 
however,  operate  as  planned  after  the  boost  phase.  We  ascribe  this  failure 
to  a  malfunction  of  the  high  voltage  supply.  The  high  voltage  supplies  of 
units  ASE-P-32-2,  ASE-P-32-3,  and  ASE-P-32-5  were  changed  to  a  new 
supply  of  the  type  used  on  the  P-11  instrumentation. 

Unit  ASE-P-32-5  is  presently  on  board  of  an  orbiting  spacecraft 
launched  on  9  Mev  1963  from  VAFB  (PMR).  information  has  yet  been 
received  by  us  on  fhe  success  of  the  flight  and  no  data  have  yet  been 
made  available. 

Units  ASE-P-32-2  and  ASE-P-32-3  have  not  yet  been  assigned  a 
vehicle.  These  units  are  completely  tested  and  calibrated  and  could  be 
made  ready  for  flight  in  a  matter  of  hours. 

3.1.8  Results 

Proton-Alpha  unit  No.  ASE-P32-1  was  flown  in  an 
engineerinq  test  on  Atlas  Pod  08,  19  December  1961,  from  Cape  Canaveral. 
The  trajectory  had  an  apogee  of  1280  km,  which  took  the  Poa  into  the  inner 
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radiation  belt.  (Figure  26).  Since  it  was  primarily  an  engineering  test 
rather  than  a  scientific  data  gathering  flight,  and  since  there  was  uncertainty 
concerning  the  expected  particle  intensity  the  instrument  would  encounter, 
the  magnetic  electron  broom  was  removed,  the  field  of  view  was  increased 
and  the  energy  thresholds  were  lowered.  This  resulted  in  the  detector 
being  sensitive  to  electrons  as  well  as  protons  and  alphas. 

Simple  considerations  using  previously  known  energy  and  voltage  thres- 
threshold  settings,  the  newly  set  voltage  threshold  settinos,  and  voltage- 
energy  linearity,  leaa  to  the  following  new  ranges  of  sensitivy  for  the 
various  data  channels: 


Channel 


a 


a 


1 


2 
HP, 

i 

HP 


1 


Protons 

3.  7  >  E  >  .35  Mev 
3.  7  >  E  >  .65  Mev 
3.  7  >  E  >  1.9  Mev 
3.  7  >  E  >  2.  82  Mev 
3.  7  >  E  >  3.  15  Mev 
E  >  3.  4  Mev 
E  >  5.  3  Mev 


Electrons 
insensitive 
Insensitive 
Insensitive 
Insensitive 
Insensitive 
E  >  0.  35  Mev 
E  >  0.  6  Mev 


A  complete  calibration  was  done  for  protons.  No  calibration  v/as  done 

for  electrons.  The  detector  had  a  9°  half-angle  field  of  view  and  an  area  of 

2  2 
-  t.  25  cm  giving  a  geometric  factor  of  '  0.  08  cm  -ster.  A  sample  of  the 

type  of  data  retiiev  d  is  shov/n  in  Figures  2  7,  28,  and  29. 


Figure  27  shows  the  variation  of  the  counting  rate  In  channel  P^ 
Sensitive  to  protons  with  energy  greater  than  .  35  Mev  with  respect  to 
cost!  where  i9  is  the  angle  between  th^  axis  of  the  field  of  view  of  the 
detector  and  the  B  field.  A  similar  plot  for  protons  with  energy  greater  than 
.  65  Mev  is  shown  in  Figure  28. 


Figure-  29  shows  the  variation  of  the  rate  in  one  of  the  fast  plastic 

ors  HP,*  Si;  ce  ihe  threshold  for  coincidence 


scintillator  (Pilot  B)  disc'.. 
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TRAJECTORY  OF  POO  IN  B-L  COORDINATES 

ATLAS  POD  08  —  19  DEC.  1961 
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Figure  29 
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in  the  CsKTl)  had  been  lowered  to  accept  pulses  for  minimum  ionising 
pai tides,  this  rate  can  represent  the  rate  either  of  protons  with  energy 
greater  than  3.  36  Mev  or  of  elect; r^ns  of  energy  greater  than  .  35  Mev. 

The  second  possibility  iz  the  one  we  have  accepted  to  explain  the 
inconsistency  between  the  extrapolated  energy  spectrum  at  3  Mev  that 
can  be  calculated  by  use  of  the  low  energy  measurements  (which  fit  a 
power  spectrum  of  the  form  E  with  y  between  approximately  2.  5  and 
3.)  and  the  measured  intensities  in  HP2. 

We  conclude  that  we  have  measured  in  HP2  a  mixture  of  electrons 
and  protons  with  electrons  predominating.  We  believe  that  channels  Pj, 
P^,  cr^,  and  are  measuring  protons  with  perhaps  a  small  contam.inaticn 
of  electrons  (due  to  straggling  or  pile-up). 

The  rai-e  observed  in  HPl  is  consistent  with  it  being  sensitive  to 
protons  of  energy  greater  than  about  5.  3  Mev. 

The  counting  rate  shows  a  sharp  rise  at  an  altitude  of  -1050  km, 

B-  .29,  L  -  1.  48;  reaches  a  maximum  at  apogee,  1280  km,  B  -  .  16  and 
L  -  1.  23;  it  decreases  sharply  a-ound  700  km,  B  -  .  20,  L  ~  1.  34.  The 
interpretation  is  that  the  pod  entered  the  inner  radiation  belt  at  1050  km, 
passed  through  apogee  inside  the  belt  and  then  left  at  700  km.  These 
"edges “  when  plotted  on  a  B-L  map  correspond  to  the  general  outline  of 
the  known  natural  belts. 

Below  is  a  table  showing  the  typical  counting  rates  at  apogee.  The 
maximum  and  minimum  refer  to  the  maxima  and  minima  of  the  signal  as  it 
is  modulated  due  to  tumbling  and  to  changes  in  B,  L  with  time;  units  are 
counts  F>er  second. 
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Channel 

Max. 

Min, 

P 

1 

700 

35 

1. 

2, 

200 

1 

~  4 

~  0 

i 

~  3 

~  0 

*^2 

~  0 

~  0 

HP„ 

7000 

—  4 

IIP 

■lit.., 

5 

~  0 

Ail  of  these  Tneasurerrier.te  vvrre  obtained  from  the  count  rate  meters. 

Two  of  tiie  scalers  appear  to  have  operated  throughout  the  flight,  a-  ^  and 

2b 

KP  ;  HP  remained  constant  throughout.,  and  c..  counted  for  the  first 
2d  lo  ■»  . 

350  seconds  and  tfien  remained  constant  fci  the  rest  of  the  flight.  These 
appear  to  be  the  only  two  possible  instrumental  failures  in  the  entire  system. 

Preliminary  ccwrelations  of  counting  rate  with  the  magnetometers 
indicate  that  it  may  be  possible  to  obtain  angular  distributions  of  elections 
and  protons.  The  turubiing  occurred  in  such  a  manner  that  not  all  angles 
v.’ere  explored;  thus  it  is  impossible  to  get  complete  a.ngular  distributions, 
but  it  may  oe  possible  to  get  enough  to  be  of  significant  interest. 

We  consider  this  flight  to  have  been  successful  beyond  expectation 
botii  from  an  engineering  and  scientific  point  of  view. 

A  complete  study  of  the  angular  distributions  v/as  not  done  in  the 
first  analysis  for  we  were  mainly  interested  in  whether  or  not  the  instrument 
operated  correctly  under  flight  conditions  —  it  did,  except  for  the  two  apparent 
scaler  failures. 
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APPENDIX  C 


1.  0  BETA.  BETA.  GAMMA  INSTRUMENTATION  (CRM  11) 

This  instrumentation  v»as  dengned  to  detect  y  -rays  and  ^-rays.  As 
shovvn  in  Figure  30.  the  instrument  package  contains  two  0-ray  detectors 
mounted  with  mutually  perpendicular  detection  vectors  and  one  y  •'lay  detector. 
This  geometric  configuration  will  allow  rough  measurements  of  angular 
distributions  to  be  performed  and  correlated  with  magnetic  aspect. 

Thie  0-ray  detector  consists  of  a  plastic  scintillator  which  is  covered 

by  a  thin  aluminum  shield  and  is  viewed  by  a  ruggedlzed  photomultiplier. 

The  aluminum  shield  prevents  visible  light  from  impinging  on  the  scintillator 

and  establishes  a  minimum  energy  for  S-rays  which  enter  thu  scintillator. 

2 

The  chosen  thickness  of  6.8  mg/cm  corresponds  to  the  practical  range  of 

electrons  of  70  kev  and  for  protons  of  1.  5  Mev.  The  thickness  of  the  plastic 

2 

scintillator,  approximately  .31  g/cm  ,  corresponds  to  the  practical  range  of 
0.  8  Mev  electrons.  The  crystal  was  chosen  so  thin  to  minimize  its  sensi¬ 
tivity  to  Y-radiation.  T’e  detector  is  exposed  to  the  environment  through  a 
hole  in  the  instrument  package.  The  field  of  view  of  this  sensor  is 
approximately  1  steradian  (Figure  31).  The  detector  is  150  times  more  sensi- 
tive  to  electrons  than  to  y  -rays.  The  range  of  sensitivity  of  this  detector  is 
from  approximately  10  to  10  electrons/cm^sec.  Figure  32  shows  the  3-sensor 
unit,  and  Figure  33,  the  schematic  diagram  of  the  beta-gamma  sensor. 

The  pulse  rate  and  the  total  cuaent  from  the  detector  are  measured. 

A  two  channel  pulse  height  analyzer  is  used  to  furnish  crude  spectral  informa¬ 
tion  for  pulse  rates  from  10  to  10^  pulses  per  secona.  The  maximum  puls  j 
rate  is  determined  by  the  fixed  dead  time  of  the  counter  of  lOp  sec.  The  pulse 
rates  measure  the  number  of  pulses  of  energy  greater  than  the  threshold  energy. 
Figure  34 shows  the  beta-gamma  logic  unit,  and  Figure  35, the  schematic 
diagram  of  this  unit. 
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EXPLODED  VIEW  OF  RADIATION  EXPERIMENT 

MODEL  II  (DISCOVERER) 
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SCHEMATIC  DIAGRAM  OF  BETA-GAMMA  SENSOR 
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DIAGRAM  OF  BETA -6 
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The  current  monitor  measures  the  rate  of  energy  deposition  in  the  crystal. 

It  is  capable  of  measuring  currents  equivalent  to  energy  deposition  rates  from 
approximately  10^  Mev/sec  to  approximately  10^  Mev/sec. 

The  y  -ray  detector  consists  of  a  CsI(Tl)  crystal  scintillator  of  approximately 
4  g  weight  covered  by  a  .  3  g/cm  A1  shield  and  viewed  by  a  ruggedized  photomulti¬ 
plier.  This  detector  is  embedded  in  a  foaming  compound  which  is  used  to  support 
the  various  components  in  the  instrument  package  (Figure  36).  The  detector  is. 

therefore,  shielded  from  the  environment  by  composite  materials  roughly  equ-va- 

2 

lent  to  at  least  2  g/cm  Al.  This  scintillator  is  consequently  relatively  insensi¬ 
tive  to  electrons  of  energy  lower  than  4-5  Mev;  it  is,  however,  sensitive  to  the 
bremsstrahlung  produced  by  these  electrons  in  the  instrumentation  case. 

Therefore,  the  response  of  this  detector  can  be  used  to  distinguish  between 
incident  photons  and  electrons.  The  detector  is  sensitive  to  y -rays  of  energy 

greater  than  about  100  kev  over  4it  steradians.  The  range  of  sensitivity  of  this 

9  2 

detector  is  from  approximately  10  to  10  y  -rays/cm  sec.  This  dynamic  range  is 
obtained  as  in  the  case  of  the  electron  detectors  by  use  of  pulse  and  current 
measurements.  Figure  37  shows  the  complete  experiment,  and  Figure  38  shows  the 
various  beta  and  gamma  sensor  logic  units  inside  the  case.  Figure  39  is  the 
block  diagram  of  the  circuit,  and  Figure  4C  is  a  diagram  of  the  wiring  harness. 

1 .  1  Calibration  Procedure  for  Debris  Decay  Radiation  Instrumentation 

1.1.1  Calibration  of  Counting  Rate  Meters 

A  variable  frequency,  variable  height  pulser  producing 
pulses  of  1  microsecond  duration  was  used  to  establish  the  discriminator 
settings  of  the  upper  level  (UL)  and  lower  level  (LlJ  discriminators.  These 
setting.*?  were  adjusted  for  all  Instruments  to  be  2.2  and  1. 1  volts.  Sub¬ 
sequently  using  a  2.  5  volt  amplitude,  1  microsecond  long  pulse,  a  curve  of 
voltage  output  of  the  count  rate  meter  versus  pulse  rate  was  obtained. 
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GEOMETRY  OF  GAMMA  RAY  DETECTOR 
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Figure  37 
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BLOCK  DIAGRAM  OF  RADIATION  EXPERIMENT 

MODEL  II  (DISCOVERER) 
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The  pulse  rate  was  obtained  by  measuring  the  distance  between  pulses  on 
a  calibrated  oscilloscope  time  base.  The  method  was  identical  for  electron 
and  Y  -ray  sensor  logics. 

1. 1.  2  Determination  of  Threshold  Energies 
137 

A  Cs  source  was  used  to  determine  the  threshold 
energies  of  the  two  counting  rate  channels  of  the  electron  sensor.  A 
pulse  height  spectrum  of  the  pulses  from  the  scintillator  was  obtained  to 
determine  the  voltage  corresponding  to  the  660  kev  internal  conversion 
electron.  The  gain  of  the  photomultiplier  was  adjusted  until  0,  5  Mev 
energy  deposition  corresponded  to  a  pulse  at  the  lower  level  discriminator 
input  of  1.  1  volts.  The  upper  level  discriminator  setting  therefcffe  corres¬ 
ponded  to  approximately  1.  0  Mev  deposited.  The  measurement  was  performed 
with  a  known  source  at  a  known  distance  from  the  detector.  A  standard 
logic  was  used  with  a  standard  detector.  The  count  rate  meter  outputs  of  the 
upper  and  lower  level  discriminators  were  recorded. 

Subsequently,  for  all  electron  sensors,  the  same  source  and  geometry 
were  used,  and  the  photomultiplier  gain  was  set  by  simply  luplicatlng  the 
measured  rates  of  the  standard  unit. 

For  the  y  "ray  detector  a  similar  technique,  utilizing  the  0.  51  Mev 

22 

and  1.  28  Mev  y  -rays  from  a  Na  source,  was  employed, 

1.  1.  3  Calibration  of  Logarithmic  Current  Meters 

The  y  -ray  detector  rate  of  energy  deposition  channel 
60 

was  calibrated  by  use  of  a  Co  source  of  10  milllcurles.  A  standard 
detector  was  placed  at  various  known  distances  with  respect  to  the  source, 
and  the  various  anode  currents  were  measured.  With  the  same  detector 

at  the  same  gain,  the  charge  deposited  at  the  anode  by  a  single 
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0.  5!  Mev  y  -ray  from  the  Na  source  was  measured.  The  charge  Oj 
deposited  at  the  anode  by  1  Mev  was  then  calculated.  The  rate  of  energy 
deposition  in  the  standard  detector  at  each  position  was  derived  by 
the  equation 

V  -1 

^1  “  'o 

I 

In  this  manner  a  tabulation  of  energy  deposited  in  the  standard  y  -ray 
detector  as  a  function  of  distance  was  obtained.  All  y  -ray  detectors  were 
then  calibrated  by  measuring  the  current  output  at  each  of  the  points  at 
which  the  rate  of  energy  deposition  was  previously  determined.  A  graph 
of  voltage  output  versus  Mev  deposited  per  second  for  each  y  -ray  aetec- 
tor  was  thus  obtained. 

The  calibration  of  an  electron  detector  rate  of  energy  deposition 

channel  proceeded  in  the  following  manner.  The  standard  y  -ray  sensor 

was  used  with  the  electron  logic  to  be  calibrated.  The  logarithmic  current 

meter  signal  output  at  each  of  the  standard  positions  with  respect  to  tlie 

Co^^  source  was  then  recorded.  Since  all  electron  and  y  -ray  sensors 

had  gains  adjusted  so  that  a  given  energy  deposited  by  v  single  event 

In  the  crystal  yielded  the  same  pulse  height  into  the  discriminator,  it 

can  be  shown  that  the  rates  of  equivalent  energy  deposition  R  ,  in  the 

60 

electron  sensor  crystal  at  the  standard  positions  of  the  Co  source  were 
given  by 


c 
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where  is  a  correction,  experimentally  determined,  due  to  the  different 
characterlatic  decay  times  of  the  CsKTl)  and  the  Pilot  B  •cintillatore  and 
to  the  properties  of  the  electronic  circuitry.  This  method  was  used  because 
it  was  not  practical  at  the  time  to  use  a  10  milUcurle  electron  source  and 
because  of  the  relative  low  sensitivity  of  the  electron  detector  to  \  -rays. 

A  calibration  curve  of  logarithmic  current  meter  output  versus  energy 
deposited  was  thus  obtained. 

1.1.4  Calibration  for  Very  Large  Fluxes 

At  fluxes  in  excess  of  10^  particles  per  second  or 
10^  Mev  deposited  per  second,  both  the  counting  rate  meters  and  the  rate 
of  energy  deposition  channel  (logarithmic  current  meter)  became  saturated. 
In  such  situations,  the  photomultiplier  used  all  of  the  current  that  the 
high  voltage  supply  could  furnish.  The  high  voltage  therefore  dropped  as 
the  current  tried  to  Increase  above  the  power  supply's  rated  rnaximum  in  a 
manner  depending  on  the  regulating  circuitry  of  the  high  -'oltage  supply. 

One  could,  therefore,  by  monitoring  the  high  voltage,  measure  the  rate  of 
energy  deposition,  thereby  increasing  the  dynamic  range  of  the  instrument 

9 

by  iibout  3  orders  of  magnitude  up  to  10  Mev  deposited  per  second.  This 

calibration  proceeded  as  follows.  Since  a  1  curie  source  was  not  available 

(which  would  have  been  necessary  to  produce  the  desired  fluxes),  the 

amounts  of  incident  visible  light  equivalent  to  various  high  rates  of  energy 

deposition  in  the  crystal  were  established.  The  intensity  of  a  visible 

3 

light  source,  covered  with  neutral  density  of  10  attenuation,  was  varied 
until  the  reading  of  the  logarithmic  current  meter  at  the  highest  energy 
deposition  not  yet  saturating  was  reproduced.  By  removing  filters  of 
Jcnown  attenuation,  intensities  corresponding  to  known  higher  rates  of 
energy  deposition  could  thus  be  obtained.  This  calibration  was  performed 
on  the  standard  unit.  For  all  other  units  it  wae  assumed  that  the  high 
voltage  supply  and  monitor  would  yield  identical  results.  This  assumption 
was  checked  in  a  few  samples  and  found  to  hold. 
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1.1,5  Use  of  Standard  Units 

All  of  the  calibrations  were  accomplished  by  comparison 
of  the  units  to  be  Cdilbrated  with  the  standards.  This  constitutes,  of  course, 
a  relative  calibration.  Absolute  calibrations  were  performed  on  the  standard 
units  above.  These  units  have  been  retained  so  th^t  if  necessary  the 
absolute  calibration  procedure  can  be  refined  and  utilized  for  these  data  for 
all  units.  ' 

The  technique  of  measurement  consists  of  the  use  of  scintillation 
counters  as  detectors  of  the  radiation  field.  A  scintillation  crystal  .esponds 
to  the  incoming  radiation  .  y  producing  light  in  quantities  proportional  to  the 
amount  of  energy  lost  by  ionizing  radiation  in  the  crystal.  This  light  is 
detected  by  the  photosensitive  surface  of  a  photomultiplier  from  which  a 
number  of  electrons  proportional  to  the  Incident  light  flux  is  produced.  The 
electrons  are  then  accelerated  from  dynode  to  dynode  of  the  multiplier  struc¬ 
ture  v/ith  a  total  number  gain  at  the  collecting  electrode  (anode)  of  approxi¬ 
mately  one  million.  Single  ionizing  particles  transversing  the  crystal  will 
produce  a  pulse  of  electrons  at  the  anode  whose  duration  is  a  function  of 

the  crystal  properties.  In  CsI(Tl)  this  duration  is  approximately  10 

*8 

seconds,  Pilot  B,  approximately  10  seconds.  Particle  fluxes  impinging 

on  tlie  scintillator  can  be  detected  by  counting  the  individual  pulses  as 

6  8 

long  es  their  rate  is  less  than  10  in  CsI(Tl)  or  10  in  Pilot  B.  If  fluxes 
larger  than  these  occur,  then  an  average  current  Is  produced  at  the  anode 
which  is  a  direct  measure  of  the  rate  of  energy  deposition  in  the  crystal. 
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To  utilize  fully  the  information  furnished  by  the  detector,  the  rate 
of  individual  pulses  as  well  as  the  average  rate  of  energy  deposition  in 
the  crystal  was  measured.  The  maximum  rate  at  which  single  pulses  could 
be  counted  was  dc‘«‘rmined  by  the  electronics  and  was  of  the  order  of  10^ 
per  second.  At  the  rates  which  permitted  pulse  techniques  to  be  used, 
pulse  height  analysis  in  two  channels  of  different  thresholds  was  per¬ 
formed  to  obtain  rough  information  on  the  spectrum  of  the  Impinging 
radiation.  For  incident  fluxes  greater  than  10^,  measurementc  were  made 
of  the  rate  of  energy  dspositlon.  These  combined  techniques  permitted 
measurements  over  a  range  from  about  10  Mev  deposited  per  second  to 

9 

10  Mev  deposited  per  second. 

In  order  to  carry  out  the  proposed  experiment  to  measure  the  flux  of 
electrons  and  y  -rays  from  radioactive  debris  decay,  two  crystals  were 
used  --  one  of  CsI(Tl)  and  the  second  of  Pilot  B  plastic.  The  reasons  for 
this  choice  were  as  follows:  since  scintillator  crystals  respond  tc  ionizing 
particles,  rheir  response  to  n  -radiation  occurs  primarily  by  creatior  in 
the  crystal  of  Compton  electrons  and  detection  of  these  electrons.  In 
order  to  have  an  efficient  detector,  one  m.ust  therefore  furnish  enough 
mass  for  the  y  -ray  to  Interact  and  a  material  of  high  atomic  number  7. 
to  provide  a  high  Compton  cross-section.  On  the  other  hand,  )f  a  low 
Z,  small  mass  cr>'s.c;'.  is  chosen,  lis  efficiency  to  \  -radiation  is  small. 

The  chosen  C'I(Tl)  and  Pilot  B  scintillators  had  messes  of  approxi¬ 
mately  4.  u  and  0.  31  grams  and  the  ratio  of  the  effective  atomic  numbers 
was  of  the  order  of  10.  Therefore,  the  CsKTl)  was  approximately  two 
orders  of  magnitude  morr  efficient  in  detecting  y  -lays  than  the  Pilot  B 
scintillator.  Both  cr/it^uls  were  approximately  equally  efficient  fvjr 
detection  of  eiecUens.  By  shielding  the  CsKTl)  crystal  with  about 


/  2 

2.0  g/cm  of  aluminum  and  composite  materials,  the  sensitivity  of  the 

CsI{Tl)  detector  to  electrons  was  reduced  by  orders  of  magnitude  for 

electrons  with  energy  less  than  about  4  Mev.  The  shielding  in  front  of 

the  Pilot  B  scintillator  was  made  as  thin  as  possible  and  was  of  0. 165 
/  2 

g/cm  of  beryllium.  This  corresponded  to  an  electron  energy  threshold 

of  0.  5  Mev.  In  Figure  41  the  calculated  fraction  of  Incident  energy 

deposited  in  the  scintillator  as  a  function  of  the  energy  of  the  incident 

radiation  Is  shown  for  the  electron  detector  for  y  -rays  and  electrons. 

The  calculation  was  carried  out  taking  into  account  the  transmission  of 

the  beryllium  filter  and  the  stopping  power  of  a  thickness  of  carbon 

equivalent  to  the  thickness  of  the  Pilot  B  scintillator.  In  Figure  42  a 

similar  plot  Is  shown  for  the  >  -ray  detector.  An  absorber  thickness  of 
2 

2  g/cm  of  aluminum  was  assumed  and  the  stopping  power  of  the  CsI(Ti) 

2 

crystal  of  about  4.5  g/cm  was  taken  into  account..  Figure  43  shows,  for 
a  given  incident  gamma  ray  fiur  the  ratio  of  the  e.nergy  deposited  in  the 
gamma  ray  detector  to  the  energy  deposited  in  ihe  electron  detector  as  a 
function  of  the  gamma  photon  energy 
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PERCENTAGE  OF  FLUX  ABSORBED 


PERCENTAGE  OF  IMPINGING  FLUX 
ABSORBED  IN  ELECTRON  DETECTOR 

VS. 

ENERGY  OF  INCIDENT  RADIATION 


0.1  1.0  to 

INOOENT  ENERGY  (Mmv) 


Figure  41 
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PERCENTAGE  OF  FLUX  ABSORBED 


PERCENTAGE  OF  IMPINGING  FLUX 
ABSORBED  IN  GAMMA  RAY  DETECTOR 

VS 

ENERGY  OF  INCIDENT  RADIATION 


Figure  42 
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RATIO  OF  ENERGY  DEPOSITED  BY  GAMMA  RAYS 
IN  GAMMA  RAY  TO  ELECTRON  DETECTOR 

VS 

ENERGY  OF  INCIDENT  RADIATION 


APPENDIX  D 


DIGn’AL  COMPUTER  PROCESSING 


1.0  INTRODUCTION 

The  telemetered  data  from  the  satellite  was  digitized  md  initially 
processed  by  loclcheed.  This  processing  produced  digital  magnetic  tapes, 
in  a  format  compatible  for  reading  by  the  FORTRAN  II  input -output  system 
on  the  IBM  7090.  These  tapes  contained  voltage  readings  telemetered 
from  the  satellite,  with  universal  times  of  original  instrument  reading  and 
mergeo  ephemeris  information.  During  the  course  of  the  satellite's  life¬ 
time,  Lockheed  changed  the  format  of  these  tapes.  The  two  formats  are 
referred  to  in  this  document  as  "Fonnat  1"  and  "Format  2".  Format  2 
gives  information  arranged  very  differently  from  Format  1.  Therefore, 
different  computer  programs  had  to  be  rewritten  for  the  two  formats. 

To  aid  in  the  reduction  of  the  data,  the  following  computer  programs 
were  written.  These  programs  all  run  on  the  IBM  7090,  under  control  of 
the  FORTRAN  Monitor  System, 

1.  Tape  Inventory  Program,  Format  1 

2.  Tape  Listing  Program,  Format  1 

3.  Count  VI.  B  and  L  Program.  Format  1 

4.  Time  Correction  Program,  Format  1 

5.  Alpha  VS,  N  (Pitch  Angle  vs.  CounO  Program,  Format  1 

6.  Alpha  vs.  N  (Pitch  Angle  vs.  Count)  Pregram,  Format  2 


The  first  program  was  designed  to  aid  in  bokkeeping  during  the  process¬ 
ing  of  the  data.  It  produced  a  summary  of  the  information  on  the  Lockheed 
tapes.  The  second  program  produced  a  detailed  listing  of  selected  informa¬ 
tion  from  the  Lockheed  tapes.  The  third  program  mapped  count  rate  informa¬ 
tion  in  geomagnetic  B-L  space.  The  fourth  program  applied  a  time  correction, 
which  had  been  computed  manually,  to  the  information  on  the  Lockheed  tape. 
The  last  two  programs  tabulated  voltage  readings  with  their  associated  count 
rate  values  and  pitch  angles. 

Since  most  of  the  information  received  wcs  in  Format  1 ,  most  f  the 
programs  ’.vere  net  rewritten  for  Format  2. 
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2.0  LOCKHEED  TAPE  FORMATS 


The  following  information  applies  to  both  formats: 

1.  Real  time  channels  are  8,  11.  12. 

2.  Tape  recorded  channels  are  18  and  98  (F) 

3.  Telemetry  Link  1  is  used  for  launch  only. 

4.  Telemetry  Link  2  contains  all  experiment  aata  except  SAPUT 

5.  Telemetry  Link  3  is  SAPUT  only. 

6.  The  station  code  is: 

01  VAFB  23  Ascention  Isia.-jd 

07  Kodiax,  Alaska  30  New  Boston,  N.  H. 

09  Hawaii 


2.  1  Format  1 

1.  The  records  and  files  on  the  tape  are  arranged  in  the 
usual  FORTRAN  2  manner,  256  words  per  record,  and 
as  many  records  as  necessary  in  a  file. 

2.  Each  file  contains  any  nurrher  of  records.  The  first 
record  contained  2  control  words  followed  by  3  data 
words  (the  rest  of  the  256  words  are  not  used  in  the 
record).  The  first  word  is  a  FORTRAN  2  control  word. 
The  second  word  tells  how  many  data  words  are  follow¬ 
ing  in  the  record.  The  next  three  give  tiie  vehicle 
number,  channel  numoer.  telemetry  line  number,  cow- 
mutator  position,  orbit  number,  station  number,  ani 
computer  run  number.  For  example,  this  Indicitive 
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might  be:  140108020020010091.  1401  vehicle  number; 

08  channel  number;  02  telemetry  link  number;  002  commutator 
position;  001  orbit  number;  09  station  number;  and  1 
computer  nin  number. 

3.  The  second  and  succeeding  records  in  the  file  will  con¬ 
tain  the  following:  The  first  word  in  the  record  is  a  FOR¬ 
TRAN  2  control  word,  and  the  second  word  tells  how  many 
words  of  data  follow  in  the  record.  Then  follow  36  {or 
fewer)  sets  of  words  arranged  as  TIME,  GG  LATITUDE. 
LONGITUDE,  ALTITUDE  (KM),  B.  L  These  are  arranged  in 
groups  of  seven.  Each  group  gives:  time,  voltage,  read¬ 
out,  geocentric  latitude  in  degrees,  terrestrial  east  longitude 
in  degrees,  altitude  in  kilometers,  and  geomagnetic  coor¬ 
dinates  B  (gauss),  an  L  (earth  radii). 

4.  Each  file  contains  all  the  data  from  one  commutator  point 
acquired  from  one  pass.  In  general  the  first  two  words  in 
each  record  are  a  FORTRAN  control  word  and  the  number  of 
data  words  following  in  the  record.  Time  is  the  common 
base  in  relating  one  file  to  other  files  and  is  given  to  two 
decimal  places,  in  seconds  of  Universal  Time  since 

1  October  1962. 

2.  2  Format  2 

The  tape  contains  any  number  of  flies.  Each  file  contains: 

1.  Flight  Indicative  in  the  first  record. 

2.  Data  indicative  in  the  second  record. 

3.  Actual  data  values  from  the  third  up  to  the  last  record 
prior  to  an  end  of  file. 
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There  Is  an  end  of  file  after  each  file  of  data  and  throe  (3)  end  of  files 
after  the  last  file  of  data  on  the  tape.  The  following  read  statement  may 
be  used  to  read  all  records: 

READ  TAPE  IT,  IRM,  MG.  NSET,  [  (DATA(I ,  J),  I  »  1,  MG),  J  «  1,  NSET] 

where 

IT  *  Input  Tape  Unit  Number 

IRM  *  which  type  of  record  this  is: 

If  IRM  •  0.  the  record  contains  three  (3)  words  Flight  Indicatives. 
If  IRM  -  1,  the  record  contains  Data  Indicatives.  This  shows 
the  arrangement  of  the  data  in  all  the  following 
data  records. 

If  IRM  -  9,  thfj  record  contains  the  actual  data  values. 

MG  =  Maxi  mum  nur  oer  of  the  first  subscript  I  in  Data  (I,  J)  of 

the  read  statement.  This  value  is  always  equal  to  the  number 
of  merged  quantities  in  the  record. 

NSET  =  Maximum  number  of  the  second  subscript  J  in  Data  C,  J). 

Normal  data  records  except  the  last  partial  record  have 
NSET  =  60.  Therefore,  the  dimension  requires  Data  (60,60). 

MG  and  NSET  are  used  by  the  Read  Tape  statement  to  read  one 
complete  record  from  the  input  tape. 

DATA  (I.J)  =  Locations,  where  the  actual  Indicative  Information  or  data 
values  are  stored. 

The  format  can  best  be  described  by  an  example: 

Suppose  a  file  of  data  to  be  read  in,  from  the  input  tape,  contains  the 


following  data: 


A.  night  1401 

Channel  8 

Link  2 

Orbit  30 

Station  30 

Run  Number  1 


B.  Pin  #3,  #4.  #11.  #38,  and  #41  data  from  Channel  8, 
Link  2  commutator  are  written  together  with  ephemerls 
data  on  each  data  record. 

After  reading  each  record,  test  IRM  value. 

1.  If  IRM  =  0,  the  Flight  Indicative  gives  BCD  information 


as  follows; 

Flight  Ch 

Link  Pin 

Orbit 

Sta.  Run 

1401  08 

- - - ( 

020000 
- - - * 

0  3  0 
^ - 

3  0  1 

Data  (1.1)  Data(2,  1)  Data  (3.  1) 


This  record  contains  all  the  information  in  Item  A  above. 

MG  s  3, 

NSET  = 

1  in  this  case. 

If  IRM  = 

1,  Data  Indicative  gives  BCD  information  as  follows 

DATA  (1 . 

1)  = 

TIME** 

DATA  (2, 

1)  = 

2***** 

DATA  (3. 

1) 

4***** 

DATA  (4. 

i)  = 

11**** 

DATA  (5. 

1)  = 

36**** 

DATA  (6. 

1) 

4  ^  **** 

DATA  (7. 

1)  » 

GC*LAT 

DATA  (8, 

1)  « 

LONG.  * 

DATA  (9,  1)  • 

ALT.  ** 

100 


DATA  (10.  1) 
DATA  (11.  1) 


•  *][^*** 


This  record  should  give  all  the  Information  on  the  order  of 
all  merged  data  In  Item  B.  MG  *  11.  NSET  ■  1. 

If  IRM  =  9,  this  record  contains  the  actual  data  which  looks 
like:  (all  data  between  the  dotted  lines  are  omitted  for 


convenience) 


DATA  (1.1) 
DATA  (2,  1) 
DATA  (3.  1) 
Data  (4.  1) 


13209,05 
1.46 
4.35 
18.  40 


Time  in  seconds. 
''In  #2  data  value. 
Pin  #4  data  vulue. 
Pin  #1 1  data  value. 


Merged 
Data  Set 


DATA  (11.  1) 


9.34 


Ephemerls  L  value 


Record  I 
of 

Data  DATA  (1.2) 
DATA  (2.  2) 


13210.05 

2.08 


Time. 
Pin  #2 


Second 
Merged 
I  Date  Set 


DATA  (11.  2) 


15.03 


L  vslue 


DATA  (1.  NSET)  ~- 
(2.  NSET)  - 
(11. NSET)  - 


14093.00 
1. 68 


Time 
Pin  #2 

L  value  of  the  last 
set  in  itwi  reccrS. 


/ 

MG  'll-  -  (one  tl.ne  word)  ♦  (5  data  pin  words)  *  (5  ephemerla  words) 
NSET  «  60  --  (last  record  may  b?  partial) 


3.0  LOCKHEED  TAPE  INVENTOPY  PROGRAM  (FORMAT  1) 


This  program  products  a  printed  summary  of  the  contents  of  the  Lockheed 
tapes.  For  each  file,  the  program  prints: 

1.  Summary  of  the  indicative  record.  This  gives  the  3  data  words  in 
the  record  in  their  original  packed  BCD  form,  and  also  in  an 
unpacked  form  with  labels  to  show  the  mea.nings  of  each  quantity 
given.  If  the  indicative  record  contains  more  than  three  data 
words,  the  program  Drint.»  an  enro.  comment  and  gives  the  first 
few  words  converted  to  floating  point  decimal. 

2.  Total  number  of  words  and  records  (not  counting  the  indicative 
in  the  file. 

If  any  record  afier  the  first  in  the  file  contains  3  or  fewer  data  words, 
the  program  assumes  that  this  is  intended  to  be  the  start  of  a  new  iile,  prints 
an  error  comment,  and  treats  this  record  as  un  Indicative  record. 

The  program  can  process  up  to  3  Lockheed  tapes  simultaneously.  It 
cannot  handle  two-  reel  'es,  but  v/ili  handle  the  second  reel  as  a  separate 
Lockheed  *dpe.  This  simply  means  ‘het  the  first  file  on  the  tape,  which  wi’l 
probably  be  a  contlncatior  of  c  file  on  the  previous  tape,  will  be  treated  as 
an  independent  file. 

The  program  uses  no  sense  switches:  It  reads  one  input  card,  which 
has  data  fields  in  columns  1-6.  /-12,  and  13-18.  In  each  of  these  fields 
the  user  may  punch,  right  adjusted,  the  logical  tape  number  whe.'‘e  a  Lockheed 
tape  will  be  mounted.  The  first  blank  or  zero  field  on  the  card  will  signal 
the  end  of  t.he  list  of  tape  numbers.  Columns  19-24  must  be  blank. 


4.0  LOCKHSED  TAPE  USTING  PRCXSRAM  (FORMAT  1) 

This  program  is  designed  to  produce  a  summary  of  the  contents  of  the 
Lockheed  tape  that  Is  more  detailed  than  that  produced  by  the  tape  inventory 
program.  like  the  tape  Inventory  program,  it  uses  no  sense  switches,  and 
reads  as  input  a  single  card,  which  is  identical  to  the  card  read  by  the 
tape  inventory  program.  For  each  file,  it  produces  the  following  output; 

1.  A  summary  of  the  indicative  record.  This  summary  is  identical 
to  that  produced  by  the  tape  inventory  program, 

2.  Every  30  records  in  the  file.  It  prints  out  the  time,  ephemeris, 
and  voltage  Information  for  the  first  5  voltage  readouts  in  the 
record.  These  groups  of  t>  voltage  readouts  are  separated  by 

a  blank  line. 

Like  the  tape  inventory  program  this  program  can  detect  the  occurrence 
of  indicative  records  in  the  middle  of  data  files. 
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5.0  COUNT  VS  B  +  L  PROGRAM  (FORMAT  1) 


1  fjtfP9S,S- 

The  program  processes  data  from  the  Lockheed  satellite  tapes 
(Format  1)  to  produce  the  following  printouts: 

1 .  A  group  of  tables  showing  the  calibrated  voltage  readouts  (count 
numbers)  as  received  from  the  satellite  as  a  function  of  the 
geomagnetic  coordinates  of  (B  and  1)  the  satellite.  Each  of 
these  tables  represents  one  voltage  point  or  commutator  position, 
and  one  range  in  B  and  L.  This  Information  is  specified  in 
inout  cards.  So  is  the  calibration  tables  used  to  convert  voltage 
to  count  number.  Any  number  of  these  tables  may  be  produced 
per  computer  run,  but  only  nine  will  be  produced  during  one  pass 
through  the  Lockheed  tape. 

2.  For  each  of  the  tables  above,  an  accompanying  table  In  the 
same  format  showing  the  number  of  points  averaged  to  produce 
each  quantity  in  the  above  tables. 

(1)  and  (2)  above  are  referred  to  as  "Table  2''. 

3.  A  straight  tabulation  of  information  on  the  Lockheed  tape.  This 
lists  time,  geocentric  coordinates,  geomagnetic  coordinates,  and 
raw  voltage  readouts  for  requested  commutator  positions  during  a 
requested  time  Interval.  This  is  referred  to  as  "Table  I".  It  is 
generally  not  requested. 

4.  A  listing,  for  a  requested  Item  in  Table  2.  of  the  count  numbers 
and  associated  times  used  tc  produce  this  average. 
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5. 


Other  information  that  is  useful  in  interpreting  the  run.  This 
includes: 


a)  Tape  inventory,  a  summary  of  the  files  on  the  tape. 

b)  Monitor  performance,  summarising- the  occurrence  of 
noise  as  indicated  from  monitor  voltages. 

c)  Summary  of  handling  of  voltage  readouts  during  Table  2 
processing.  This  tells  how  many  voltage  readouts  had 
B  +  L  values  that  were  outside  the  range  of  the  B-L 
table,  and  how  many  were  rejected  due  to  monitor  failure. 

6.  When  there  are  errors  in  the  Lockheed  tape  or  in  the  input  cards, 
comments  lo  enable  diagnosis  of  the  errors. 

While  producing  Table  2,  the  program  rejects  all  points  during  times 
of  noise,  as  indicated  from  monitor  voltages  which  should  remain  steady. 
Monitor  voltage  points  and  nominal  voltages  are  indicated  in  the  input  cards. 

The  Lockheed  tape  may  be  in  one  or  two  reels. 

5.  2  Functional  Description 

The  program  is  essenti.\lly  a  straightforward  data  processing 
program,  and  the  functions  that  it  perfom.s  are  described  in  Section  4.  0 
in  the  text  of  this  report.  However,  the  user  shoula  be  aware  of  a  few  addi¬ 
tional  points; 

1 ,  Conversion  from  voltage  readout  to  count  rate  is  performed  by 
linear  interpolation  using  the  calibration  curves  supplied  as  input. 

2.  Voltage  readouts  are  considered  by  the  program  to  be  '‘noi;ie-free" 
1^  no  monitor  voltage  failure  has  been  detected  within  T  seconds 
of  the  time  of  the  voltage  readout.  The  value  of  T  is  read  as 
input. 
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3.  Each  B“L  box  In  Table  2  represents  the  average  of  count  rotes  that 
fell  within  this  B-l  interval.  To  compute  each  individual  count 
rate  value,  the  average  of  5  succ^slve  voltage  readouts  is 
used.  When,  due  to  occurrence  of  noise,  less  than  5  consecutive 
readouts  are  available,  then  these  readouts  are  not  used. 

5.  3  Operating  Instructions 
Sense  Switches; 

ON  -  If  2  binary  tapes  are  in  use; 

OFF  -  If  i  binary  tape  Is  in  use. 

Tapes: 

Tape  A-6  is  the  binary  tape  prepared  by  Lockheed  which  i.icludes 
satellite  ephemeris  in  both  geographic  and  geomagnetic  coordinates,  plus 
voltage  readouts.  If  this  tape  Is  in  two  reels,  the  second  reel  must  be 
mounted  on  unit  B-6,  and  sense  switch  2  must  be  on. 

Input  Cards: 

Tape  A-2  is  prepared  from  a  card  deck,  which  may  be  stacked  with 
other  FMS  jobs.  The  makeup  of  this  deck  is  as  follows: 

FMS  ID  curd* 

XEO  card* 

PROGRAM  BINARY  DECK 
FMS  DATA  CARD  * 

Deck  of  control  cards  read  by  the  program 
The  control  card  deck  tells  the  proyran  what  data  to  process  from  the 

*  These  are  standard  FMS  control  cards. 


1C6 


Lockheed  tape,  and  how  to  process  it.  It  contains  the  following  sets  of 
cards,  in  the  following  order: 

1.  Channel  specification  cards 

2.  Table  2  specification  cards  {There  may  be  any  number  of  sets 
of  these  after  a  set  of  chanrel  specification  caids. 

3.  Recycle  card 

The  sequence  above  Is  referred  to  as  a  Channel  Deck,  and  may  be 
repeated  any  number  of  times. 

Set  (1),  Channel  specification  cards,  contains  the  following  cards: 

1.  1  Page  Heading  card 

2.  1  Channel  Request  card 

3.  1  Taole  1  Request  card 

4.  Up  to  10  Monitor  Specification  cards.  [  The  last  Monitor 
Specification  card  is  differentiated  by  the  presence  of  data 
in  the  4th  fit'd,  which  must  be  blank  in  the  previous  Monitor 
Specification  cards,  (See  Card  Formats,  below).  ] 

Set  (2),  Table  2  specification  cards,  contains  the  following  cards: 

1.  1  Table  2  Request  card.  (This  gives  the  commutator  position 
number  and  limits  for  the  B-L  table). 

2.  A  block  of  cards  which  specify  the  Calibration  Curve 
This  block  of  cards  contains: 

a)  Calibration  Curve  Size  card 

b)  Enough  Calibration  Curve  Table  cards  to  represent  all 
the  voltages  (ordinates  of  the  curve) 
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c)  Enough  Calibration  Curve  Table  cards  to  represent  all 
the  particle  counts  (abscissas  of  the  .-’.urve). 


This  block  of  cards  may  be  omitted  altogether  if  the  following  conditions 
are  met: 

a)  This  same  calibration  curve  was  given  previously  for  this 
commutator  position. 

b)  No  tscycle  card  has  appeared  since  giving  the  calibration 
curve. 

3.  1  B-L  Box  Detail  Request  card  (This  specifies  number  of  B-L 
box  lists  to  be  printed  01*1.) 

4.  N  B-L  Box  List  Request  cards,  where  N  *  the  number  of  lists 
specified  in  the  B-L  Box  Detail  Request  card.  These  cards  must 
be  omJited  if  the  B-L  Box  Detail  Request  card  is  blank. 

Card  Formats; 

The  above  formats  of  the  above  cards  are  given  below: 


Page  Heading  Card 


Columns  1  -  72: 


Channel  Request  cart'i 
Columns  1-2: 


3  -  12: 


Anything  punclied  here  will  appear  as 
the  first  line  of  each  page  of  output. 

It  is  suggested  that  Identifying  informa¬ 
tion  from  the  tape  label  be  punched  here. 

Channel  numbers  of  data  desired  from  the 
binary  tape.  This  number  must  be  right 
adlutted  and  filled  In  with  a  leadina 
zero  if  appropriate. 

Must  lie  blank. 
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13  -  22:  Start  time  for  printed  voltage  tabula¬ 

tion  (Table  1). 

23  -  32:  Stop  time  for  printed  voltage  tabula¬ 

tion  (Table  1).  These  times  must 
be  in  the  same  units  as  used  on  the 
binary  tape,  and  must  be  punched  with 
a  decimal  point,  anywhere  in.  the  field. 

33  -  40;  Must  be  blank. 

41  -  44:  Punch  any  integer  here  to  cause  print¬ 

ing  of  messages  whenever  a  voltage 
point,  specified  in  any  of  the  Table  2 
requests  below,  has  B-L  coordinates 
outside  the  specified  P-L  table.  Leave 
blank  to  suppress  such  printing. 

Table  1  Request  Card 
Columns  1-6 
7-12 

13-18  Each  of  these  fields  may  be  punched  with 
19-24  a  right  adjusted  integer,  which  specifies 
25  -  30  a  commutator  position  whose  voltage 

31  -  36  ^  readings  are  to  be  tabulated  in  Table  1. 

37  -  42  As  many  of  these  fields  as  desired  may 

43  -  48  be  used,  but  the  first  blank  or  zero  field 

49  -  54  (sccinning  from  left  to  right)  will  signal 

5S  -  60  the  end  of  the  Table  1  requests. 

61  “  66 
67  -  72 


Monitor  Specification  Cand 


Columns  1-6: 

7  -  12: 

13  -  20: 


21  -  28. 


29  -  36: 


Ignored . 

Commutator  position  number  for  this 
monitor.  (Right  adjusted  integer) 

Nominal  voltage  reading  for  this 
monitor.  (Punch  with  decimal  point) 
Allowable  deviation  from  the  nominal 
voltage  reading.  (Punch  with  decimal 
point) . 

This  fielc  must  be  blank  for  all  but  the 
lasi  Monitor  Specification  card  in  a 
series.  In  the  last  Monitor  Specifica¬ 
tion  card  in  the  series,  this  must  be 
punched  with  a  non-zero  time  interval 
(punch  with  decimal  point)  T,  such 
that  any  data  reading  within  T  time  units 
of  a  monitor  failure  will  be  discarded. 

This  will  be  punched  as  a  negative  number 
if  it  is  desired  to  print  messages  about 
each  monitor  failure.  It  will  be  punched 
as  a  positive  number  to  suppress  such 
printing. 


Table  2  Request  Card 
Columns  1  -  i: 

5  -  10; 

11  -  16; 


Commutator  position  number  for  this 
table.  (Right  adjusted  integer) 

T4:>Wtfr  table  limit  in  B.  (Punch  with 
decimal  point) 

Upper  table  limit  in  B.  (Punch  with 
decimal  point) 


17  -  22: 

Lower  table  limit  In  L. 

decimal  point) 

(Punch  with 

23  -  28: 

Upper  table  limit  In  L. 

decimal  point) 

(Punch  with 

29  -  48: 

Name  of  point 

Calibration  Curve  Size 

Card 

Columns  1-12: 

N,  the  number  of  points  given  in  the 

calibration  Table  cards  that  follow. 
This  must  be  punched  as  a  right- 
adjusted  Integer. 


Calibration  Curve  Table  Card 


Columns  1-12 
13  -  24 
25  -  36 
37  -  48 
49  -  60 


Each  of  these  fields  may  contain  a  coor¬ 
dinate  of  one  of  the  points  on  the  calibra¬ 
tion  curve.  Numbers  should  be  punched 
with  a  decimal  point.  A  decimal  exponent 
may  be  punched,  in  the  form  E  +  xx,  and 
it  must  be  right -adjusted,  In  the  field. 
There  are  two  sets  of  these  cards  for  each 
curve.  The  first  gives  the  ordinates,  or 
voltages,  in  numerical  order,  and  the 
second  set  gives  the  corresponding 
abscissas  or  particle  counts.  Each  set 
will  rtsrt  a  new  card.  Numbers  must  be 
punched  in  successive  fields  from  left 
to  rlnht  across  the  card,  and  contloued 
onto  successive  cards  until  all  N  points 
are  represented. 


B  -  L  Box  Detail  Request  Card 

Column  1:  K^inus  sign.  (11  punch)  (May  be 

blank  If  rest  of  card  is  blank) 

^  -  12:  The  number  of  B-L  Box  List  Request 

cards  that  follow.  (Right >ad justed 
integer) . 

B  -  L  Box  list  Request  card 

Columns  1*12:  Column  number  of  B-L  box  to  be  listed 

in  detail.  (Punch  as  right-adjusted  integer). 
13-24:  Row  number  of  B-L  box  to  be  listed  in 

detail.  (Punch  as  right-adjusted  integer). 

Recycle  card 

There  are  two  possible  formats  for  this  card: 

Format  1: 

Column  4:  1 

Rest  of  Card:  Must  be  blank 

The  Format  1  Recycle  card  must  be  immediately  followed  by  another 
Channel  Deck,  unless  it  is  the  last  card  in  the  input  deck. 

Format  2: 

Column  4;  2  n 

Rest  of  Card:  Must  be  blank 

This  card  is  used  to  instruct  the  operator  to  mount  a  new  tape  on  unit  A-6 
(and  B-6  if  in  use).  It  must  always  be  followed  by  another  card,' which  may 
contain  in  columns  1-72  any  instructions  to  the  operator. 

Upon  recognising  a  Format  2  Recycle  card,  the  program  will  read  the 
next  caid,  pr.nt  its  contents  on  the  on-line  printer,  and  half.  Pressing 
start  will  cause  the  machine  to  read  another  Channel  Deck. 


Restrictions  on  Program 


Due  to  the  limited  table  sizes,  certain  restrictions  exist  and  are  given 
below: 

Restrictions  on  Monitors 

No  more  than  10  monitors  may  be  specified  in  any  one  set  of 
Channel  Specification  caids. 

No  commutator  position  higher  than  60  may  be  specified. 

Restrictions  on  Table  1 

No  more  than  12  voltage  points  may  be  listed  in  Table  1 . 

No  commutator  position  higher  than  60  may  be  specified. 

Restrictions  on  Calibration  Curves 

No  calibration  curve  may  contain  more  than  50  points. 

For  any  one  Channel  Deck  there  may  be  no  more  than  20  different 
Calibration  Curves  input.  The  total  combined  size  of  all  Calibration  Curves 
may  not  exceed  250  points. 

Restrictions  on  Table  2  Requests 

No  commutator  point  higher  than  00  may  be  specified. 

Restrictions  on  B-  L  Box  Detail  Requests 

For  any  one  Channel  Deck  there  may  be  no  more  than  20  B  -  L 


Box  List  Request  cards. 


6.0  T^CKHEED  TAPE  TIME  CORRECTION  PROGRAM  (FORMAT  l) 


This  program  is  designed  to  apply  a  universal  time  correction  to  selected 
data  from  the  Lockheed  tapes.  It  does  not  compute  the  time  correction.  That 
must  be  done  manually,  nd  given  to  this  program  as  input. 

6. 1  Operating  In*  ructions 

1.  Data  Deck 

a)  Time  correction  card,  columns  1  -  6„  +  XXXX,  where  XXXX 
is  an  integral  number  of  seconds. 

+  means  that  the  voltage  should  occij  XXXX  seconds  later  down 
the  tape. 

means  that  the  time  and  ephemeris  correspond  to  a  voltage 
occurring  .XXXX  seconds  later  down  the  tape. 

b)  Channel  selection  card,  columns  1  -  6,  XX,  channel  no. 

c)  Commutator  position  selection  card,  columns  1  -  40.  (In 
any  column  N,  a  "I"  indicates  that  commutator  position  N 
is  desired  to  be  saved. 

2  Sense  Switches 

#2  OFF,  if  Input  tape  is  I  reel. 

ON.  if  input  tape  is  2  cr  3  reels. 

#3  OFF.  if  input  tape  is  1  or  2  reels. 

ON.  if  input  tape  is  3  reels, 

3.  Tajai 

FMS  Systems  Tapes,  plus  A- 7  for  Reel  1  in  input;  A-6  for 
binary  output  (SAVE);  B-6  for  Reel  2  of  Input;  B-7  fcr  Reel  3  of  Input. 

i  1-; 


7. 0  ALPHA  VS  N  (PITCH  ANGI£  VS  COUNT  RATE)  PROGRAM.  FORMAT  1 


7. 1  Proorani  Description 

This  program  tabulates  particle  counting  rate  with  corresponding 
pitch  angle.  Input  cards  specify  what  areas  of  time  and  what  experiment 
voltage  points  are  to  be  tabulated.  A  fifth  order  divided  difference  table 
technique  is  used  to  interpolate  in  calibration  curve  tables,  to  convert  from 
voltage  readout  to  count  rate.  The  calibration  curve  tables  are  read  as  input. 

Voltage  readouts  are  considered  by  the  program  to  be  "nolse-^ef "  if 
no  monitor  voltage  failure  has  been  detected  within  10  seconds  of  the  time  of 
the  voltage  readout. 

Since  the  satellite  tumbles  with  a  period  that  is  much  less  than  the 
commutator  revolution  rate,  successive  readouts  of  the  same  magnetometer 
voltage  output  should  .show  a  slow  variation.  Therefore,  the  program  smooths 
the  magnetometer  data,  to  reduce  the  effect  of  undetected  noise  in  the  magneto¬ 
meter  readouts.  Taking  each  magnetometer  axis  independently,  it  fits  each 
group  of  7  consecutive  readouts  to  a  second  order  polynomial,  using  a  least 
squares  technique. 

These  magnetometer  readouts  are  then  used  to  compute  a  unit  vector 
describing  the  direction  of  the  magnetic  field.  Cosine  pitch  angle  is  then 
computed  as  tht*  dot.  product  of  this  vector  with  another  unit  vector,  read  as 
input.  descrit4ng  the  orientation  of  the  instrument  Ir.  the  magnetometer 
coordinate  system. 
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Data  Is  printed  in  the  following  format:  TIME,  GLAT,  GLON,  ALT, 

B,  L,  BM.  ALPHA.  COSINE  ALPHA.  VOLT.  COUNTING  PATE.  The  time  is  in 
SAf^onds  and  free-running  from  1  October, 

GLAT  Geodetic  latitude  accurate  to  1  place. 

GI/DN  Geodetic  longitude  accurate  to  1  place. 

ALT  Geodetic  altitude  in  kilometers  accurate  to  1  place 

B  Theoretical  B  field  in  Gauss  accurate  to  3  places. 

L  Geomagnetic  L  coordinates,  measured  in  earth  radii 

at  the  equator,  accurate  to  3  places. 

BM  Measured  B  field  (in  Gauss)  accurate  to  3  places. 


ALPHA  Pitch  angle. 

COSINE  ALPHA  Cosine  pitch  angle. 

VOLT  Voltage  reading  of  the  specified  commutator  point, 

accurate  to  3  places. 

COUNTING  Particle  count  rate  of  the  specificKi  commutator  point. 
RATE 


7.  2  Operating  Instructions 
Sense  Switch  2  - 

ON  if  2  reel  binary  tape  on  units  A-6  and  B-6; 

OFF  if  1  reel  binary  tape  on  unit  A-6 
Input:  Tape  A-5  Is  the  binary  tape  in  Lockheed  Tape 

Format  1:  If  it  is  a  two-reel  tape,  the  second  reel  Is  mounted  on  B-6 

Tape  A- 7  is  prepared  from  card  deck  as  follows: 

FMS  ID  Card 
XEQ  Card 
Binary  Deck 
fMS  Data  Cara 
DATA  DECK 
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The  makeup  of  the  data  deck  is  as  follows: 

1 .  Identification 

2.  Experiment  number,  order  of  accuracy,  and  speciacation 

3.  Channel  request 

4  Time  intervals 

5.  Direction  cosines  for  all  four  instruments 

6.  Magnetometer  conversion  coefficients 

7.  Calibration  curves 

8.  Monitor  data 

9.  Commutator  positions  or  magnetometer  data 


A  data  deck  may  consist  of  any  number  of  sets  of  data;  each  set  contains  the 
nine  types  of  cards  specified  above,  as  long  as  the  end  specification  on  the 
last,  and  only  the  last  set,  is  not  zero. 


Formats: 


I, 


Identification  — 
Columns 
1  -  6 
7-12 
13  -  18 
19  -  24 


Description 
Flight  number 
Orbit  number 
Station  identification 
Computer  run  number 


2.  Experiment  card 
Columns 
1  -  6 

7  -  12 

13  -  IB 


PescrlpUon 

Experiment  number  (any  number  frcm  1  -  4) 
Integer,  right -adjusted. 

Order  of  polynomial  interpolation  desired, 
integer,  right-adjusted. 

End  8pecific.*tto.n,  integer,  n  ]nl-adjusted  — 
1  if  this  is  the  last  data  set;  0  otner.vise. 


3.  Channel  request  — 

ggi.tfianf  Bi8i<gilgU9n 

1  '  2  Chennel  number  of  two  dlqlt  Integer 


4.  Time  Intervals  — 

Columns  Description 

1  12  Minimum  time, 

12  ~  24  Maximum  time. 

In  the  format  XX. .  .X.X.  there  may  be  as  many  as  10  time 
intervals  requested  in  a  data  set.  The  last  request  must  be 
followed  by  a  blank  card. 


5.  Direction  cosines  — 

Columns  Description 

1-10  X 

11-20  y  direction  cosines 

21-30  z 

In  the  format  ^  XX. .  .X.  XXX  there  are  four  such  cards,  one 

for  each  instrument.  The  order  of  the  cards  must  correspond 
to  the  order  assigned  to  the  instruments  as  experiment  numbers. 


6.  Calibration  for  magnetometer  data  — 

gglMjaai  P?t9r,ljngn 

1-10  X 

11-20  y  Factor,  in  the  foimat  XXX. .  X.XXX 

21-30  z 

Here  are  two  cards;  the  first  gives  the  numbers  which  are 

subtracted  from  the  magnetometer  voltage  and  the  second 
gives  the  division  factor. 


7. 


Calibration  curves  — 


A.  Size  Card 

S^VAUUiUg 
1  >  12 
13  -  18 

19  -  24 


l>e5CiiPtlvn 

Number  of  points  in  the  calibraticm  curve 
Number  of  the  commutator  position  for 
which  the  calibration  7  irve  follows. 
Experiment  number  for  this  commutator 
position. 


B.  Curves  — 
Columns 
1  -  12 

13  -  24 


Description 

Each  field  may  contain  a  coordinate  of  a 
point  on  the  calibration  curve. 

Numbers  should  be  punched  in  the  form 
XXX.XXE+^  and  right -adjusted  in  the  field. 
There  are  two  sets  of  these  cards  for  each 
curve.  The  first  gives  the  ordinates,  or 
voltages,  in  numerical  order,  and  the 
second  gives  the  corresponding  abscissas 
or  particle  counts.  Each  set  starts  a  new 
card.  Nurr'  ->r8  must  be  punched  in  successive 
fields  from  lei’  to  right.  A  blank  card  must 
follow  the  last  cun^e  to  be  read  in. 


Monitor  data  — 

There  are  two  cards  for  each  monltoi.  The  first  card  contains  the 
commutator  position  of  The  monltcr  as  a  right- adj  usted  Integer 
in  columns  1-4. 
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The  second  card  contiiins: 


Columns  Description 

1-4  Anything  except  0 

5-16  Nominal  value  of  the  monitor  voltage. 

17-28  Allowed  deviation  of  the  voltage.  The 

first  field  is  an  integer,  the  last  two  are 
in  the  format  X. .  .X.XXX.  The  last  monitor 
request  must  be  followed  by  a  blank  card. 


9.  Magnetometer  commutator  positions  — 

Columns  Description 

1-6  X 

7-12  y  ?  commutator  position 

13-18  z 
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Restrictions 

Only  one  channel  may  be  requested  in  any  data  set.  A  maximum  of  ten 
time  intervals  may  be  requested.  Calibration  curves  may  contain  a  maximum 
of  25  points  each. 

A  maximum  of  7  monitors  may  be  used.  Up  to  10  non-monitors,  non¬ 
magnetometer  commutator  positions  may  be  requested  per  data  set. 

If  a  total  of  more  than  2,000  seconds  cf  time  is  Included  in  the  time 
intervals,  the  program  will  have  to  make  two  extra  passes  of  the  tape. 

The  first  pass  compiles  monitor  and  magnetometer  data.  The  second  pass 
correlates  this  with  requested  commutator  positions. 

If  the  time  intervals  of  twc  succeeding  data  sets  «re  exactly  the  same, 
the  monitor,  magnetometer  pass  will  be  by-passed.  Care  should  bo  taken  so 
that  this  does  not  occur  if  the  first  of  these  data  sets  required  more  than  2.  000 
seconds,  t^cuuse  the  first  1,000  seccr.ds  of  magnetometer  data  will  be  lost 
for  the  second  set. 


8.0  ALPHA  VS  (PITCH  ANGLE  VS  COUMT  RATE)  PROGRAM.  FORMAT  2 


8. 1  Program  Description 

This  program  tabulates  voltage  readouts  from  commutator  position 
5  on  channels  11  and  98  with  cosine  of  pitch  angle.  The  output  gives 
ephemeris  quantities,  including  B  and  L  from  the  ephemeris;  value  of 
B  computed  from  the  magnetometer  readouts,  cosine  pitch  angle,  and  voltage 
readouts. 

The  program  tabulates  the  following  commutator  positions  which 
represent  experimental  data:  2,  3,  4,  6,  7,  I*"  20,  21,  23,  24,  25,  27, 

28,  29  from  channels  11  and  98,  but  does  not  give  pitch  angles  for  these 
positions.  Magnetometer  coefficients  and  direction  cosine  of  instrument  are 
read  as  input. 

8.  2  Operating  Instructions 

The  prograrr;  reads  three  input  cards.  Each  of  these  cards  has 
three  numeric  fields,  as  follows:  columns  1-10,  11-20,  and  21  -  30. 

Data  must  be  punched  with  explicit  decimal  point  in  each  of  these  three 
fields. 

The  first  card  gives  the  bias  voltage  values  for  the  x,  y,  and  z  axis 
magnetometers ,  in  that  order. 

The  secorjd  card  gives  slope  values  for  three  magnetometers.  Together, 
these  cards  allow  B  component  to  be  compvited  as  a  functrcn  of  magnetometer 


voltage,  as  follows: 


B  =  (v  -  V  )/a 
o 

where  v  is  bias  and  a  is  slope. 

0 

The  third  card  gives  x,  y,  and  :z  axis  direction  cosines  for  the  direc¬ 
tion  of  the  instrument  pointing  in  the  magnetometer  coordinate  system. 
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APPENDIX  E 


A  reprint  of  a  paper  entitled, 

"Measurements  on  Trapped  Particles  Injected  by  Nuclear  Detonations" 


L.  Katz,  D,  Smart,  F.  R,  Paoiini,  R.  Giacconi  and  P.  J.  Talbot,  Jr. , 
Space  Research  IV.  Proceedings  of  the  Fourth  International  Space 
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MEASUREMENTS  ON  TRAPPED  PARTICLES  INJECTED  BV 
NUCLEAR  DETONATIONS 


L.  KATZ  and  I>.  SMART 

Air  F(trct  C<itnhruige  Research  LaborrUorirs.  lirdford.  Mom.,  USA 

and 

F.  R  PAOLINI.  R  (ilACf’ONI  and  R  J  TALBOT,  Jr. 
Amer.can  Srirnce  attd  Engtnff.ritig,  Inr.,  Cambridge,  Mass.,  1<SA  * 


Abstract:  A  polur-orbitinff  -Vir  Force  ^uitcllitc  currifnl  aloft  a  |>ayload.  shanfl  by  many 
oxfa'rimcnlcrs,  of  ra<iiation  detection  inatnunentH  for  stiidymt;  artificially  pnxfiiced 
trappoil-radiation  lielta.  Pertinent  information  on  tin*  natcilip-'g  orbit  and  capa- 
bil'ties  i«  given.  l),ita  from  four  acint illation  detectors  arc  presi'ntrxl.  The  detectors 
have  capabilities  for  diHcriminating  aniong  tyji*>»  of  |>articl<'s.  ami  for  measuring 
t)<  ;b  energy  and  angular  diHtribiitii'na  of  trap(Mxi  particl«»a.  I>ata  shouing  spatial 
d.itribiitions  of  trappeil  radiation,  both  as  omni-direotional  flux  versus  L-shell 
plots  and  us  iso. intensity  contour  plots,  are  preseute«l.  The  major  portion  of  data 
presentetl  emphasir.es  electron  pitch  angle  distributions  and  derived  parameters. 
Bx.to  ap-'dal  and  P'lniKiral  variations  of  the  pitch  angle  distributions  (meaaurod 
with  a  6'  angidar  resolution)  are  given.  Some  tentative  conclusions  are  drawn. 

PeanMir:  Ha  cii\  iiiitKe  lllir,  artiiymeitiioM  na  no-i.ipHyio  opOHTy,  <!)U.ia  ycraHOB- 
aeiiii  amiapuTypa,  CMotiTiipoimmi's-  iiccko.iiiKhmh  HCC.ieaoBaTe.iHMH.  Ohs 
iifie.i(  TaH.iHei  codoK'  ,ieT('KTn{>w  n„..;,  jeiiiiH.  npeaioaniatPiiHue  iiH  iiayHCitHH 
ticityi  <  1  itemio  coa,,  i"’iu\  pa.tiiaHHomttK’x  iiohcob.  IIphho.ihtch  napaweTpH 
opOHTU  fliyTllilKa,  a  TUK/KC  cro  TCXIIHHe 'KHC  Xap'lKTPpHCTHKH.  IlpMHOJlHTCH 
.taiiitlje  no  MeiupcM  cumhth  i.ihuhohunm  .icTCKTopaM.  JlcxcKTopu  cnocofuiw 
pa3.iH<iaTt.  xapaKTcp  'tacTHUu  it  HaMcpHTb  KaK  aiiepPHS),  TAK  H  yr.iOBoe  pac 
iipeie.ieiiiie  auxtiimemihix  xucT'.iit.  njiHiicacHU  juhhwp,  noKiiauBUKimMe  npo- 
cTpain  Tneiiiiwe  p.ecnpe.ie.icHHH  aaxnimeniiotl  paaiiauHii:  kak  rpai|)MK  noTOKa, 
iiceiiaiip  ui.ieiiioicu  no  oTiiouietiHKi  k  ofKi.ionKe  L,  thk  h  KOif.  ypij  ofoacTert, 
OI  pailimeillllJX  IHtlHlIMH  palllloH  HIlTeHClllUKH’TH.  (k  llOBHaH  MaCTb  npHIlP- 
.leiiniix  laniiijx  .siex  pacnpe.ieieiiite  a.ieKipoiion  no  iihtm  yr.iy  it  ApyrMC 
napaMctpu.  .Iukhi  h  kak  M(‘|>(  XpHm  xneiMiue,  thk  h  BpeMeHitue  HapHauHH 
pai  npe.te  irKHli  in  IIHTM  yiviy  (H.iMepemiwe  c  yr.ionijM  paJpeiucHMCM  6‘), 
i'.le.liinu  ttCKOXopMe  npe.lMapHTe.ll-lllje  UNl<o,H.l. 

I .  IntroductioD 

III  the  fall  of  !*.(ti•J,  an  .Air  Kor<x‘  satellite  gistni  the  mteniational  dfsig 
natjoii.  MliiJ  pV,  w  it',  lautn  he<l  into  orbit  The  eat<  llite  earrieti  iibstrumeiitatioci 
to  '.fudy  artitici.M  radiatioii  !«  It^  .  n-ated  l*v  the  injeetion  of  elecfmis  front 

*  ’lim  worl,  w .M  (x  rfiinnixl  miii'  i  ('l•lllrll<  l  .\K  IB  tB2rti  J'W  with  the  \if  Kor\-c 
I'niiibridge  l<<.w>Hri  li  lAtlMirstoner.  Upsl|.>rt|.  MaMsaebunrtts 

t.4fi 
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nuclear  detonations  into  the  magnetoMphere.  Under  the  scientific  direction 
of  Air  Force  Cambridge  Research  laboratories,  a  number  of  inatruments 
were  chosen  for  inclusion  in  the  payloe-d.  In  addition  to  Air  Force  Cambridge 
Research  lAlioratories,  the  folic,  t  .  ^^inizations  participated  in  the 
experiment:  American  Science  and  Engineering,  Inc.,  Jjawrence  Radiation 
Laboratory,  Johns  Hopkins  University,  the  University  of  Utah,  Air  F 
Special  Weapons  Center,  and  I^ockheed  Missiles  and  Space  Cenip^nv 
The  pertinent  parameters  for  the  satellite  orbit  »t  early  times  were  t;  d 
following:  ajiogee,  3000  n.m.;  {Krigee,  115  n.m.;  orbit  inclination,  71°; 
orbital  period,  147  min;  apsidal  rate,  O.H°  per  day.  Fig.  1  illustrates  these 


orbital  parameters.  The  satellite  pitch  {leriod  was  !27  sec,  and  the  roll  period 
was  53  sec  The  satellite  parameters  were  well  suited  to  study  the  phenomena 
ov  interest.  For  example  the  high  eceentriviiV  of  the  e.rbit  enabled  the 
acquisition  of  data  at  large  L  vah'cs  near  ihe  equatorial  plane.  Slow  tumbling 
rates,  along  with  the  eh<wen  da^a  sampling  rate  of  atiout  one  jht  second, 
jiermittetl  the  use  of  insi. nments  which  could  rneasurr-  accurately  the  angular 
distributions  of  trap|M'd  particles. 

Fig.  2  shows  the  projection  of  an  earlv  oriiit  on  a  di5M)le  tieid.  Tlie 
coordinates  are  the  usual  "nat  H  L  nM)r(iiiiates  of  Mcllwc.in  i2]. 

The  squares  ami  circles  show  region^  of  data  acipiiMtioii  m  the  1.2  and  l.H 
L  .shells,  shells  of  mtert'Mt  with  reN|H'ct  lo  artificial  radiation  ladts 

The  instrumenf.s  comprising  the  total  pA\  lo,ad  were  e  lected  to  vieid 
n*<lun<lant  n.easun'tr.ents  in  th*'  energy  range  corresjsuiding  to  fi.vsion 
electrons,  and  D*  niciiide  mtrrsiw  n<’«.»'ptHma.  angb*  detei’tors  as  well  as  ornni 
directional  deterUon  I iistrtimentN  obtained  their  |s>wer  from  a  large  capacity 
bat»ery  pack  to  avoid  degra<!»li<»n  of  the  esjierimetit  due  l<i  possible  solar 
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Kijj.  ‘J.  of  an  orbit  in  H-L  apaoo, 


coll  (iarnagf.  Data  could  lx*  acqiiiri'd  ladh  by  direct  ‘‘real-time"  read-out 
and  by  coitiniand  read  out  of  a  ta|K*  recorder  capable  of  storing  data  from 
all  ju.struniciits  for  a  time  Iouu<t  than  the  orbital  |)triod. 

Tlii-  pa|)er  will  deal  with  the  data  obtaiiu d  from  the  AFCRL-ASE  portion 
of  the  payload  only.  Data  obtained  by  other  exf)eri!nenter8  are  reported 
el.-cewhere  [s].  At  this  time,  only  a  small  jMirtion  of  the  available  data  has 
tKa-n  analyze<l.  This  is  esja-oiaily  true  of  data  referring  to  .spectral  information. 
In  A  previous  pajK-r  [1]  the  authors  pre.s<u\ted  data  illustrating  the  gross 
a.s{»<-ct.s  i  f  th(  artificial  belts.  This  pajw  (  luphasizes  the  measured  pitch 
angle  di- tributions  as  functions  of  both  time  and  magnetic  shell,  .nnd  the 
intormatio  I  derived  from  tlu  ni  Kfibrts  of  armlysi.s  have  Iktu  concentrated 
on  data  obtained  lK-twe«>n  27  OetolxT  and  30  Oetolier.  As  a  re.sult.  the  effects 
of  primarily  a  .single  nuclear  detonation  a  e  investigate!],  allxdt  fairly 
(•on;pl“ieiv. 

i’rogrimv  fur  both  Topo  and  H‘>2u  IBM  eom[Miters  have  t)e<*n  develor«*d 
ami  ns«‘d  ulienevir  [Miv'ibb-  in  data  n'dnction.  Digitized  tajH*s  with  merged 
epiu  111.  rides  n  eluding  ft  /.  e<M)nl(  lates  derived  fn»m  McIIwain’s  imslel  [2], 
and  pnm  Iml  i  anls  Iw. icing  detect. >r  ealiiirations  an-  primary  data  inputs  for 
such  program'.  Untput'  melu<l*>  a  H  L  matrix  with  deN-etor  (smuting  rates 
a*  element.,  I'rngrams  abo  e\i-t  for  yieliling  jutch  angle  distrihiitiiuis 
dinetly  wbe.i  in.igni  l< itm  ler  tlata  an*  addeil  t*i  the  input 

2  Instrumentation 

rill-  In  f  'lint  I’i.it  |i  ift  detiefiU"  einnprisiug  the  AM'Bl.  .-VSE  |s>rtli'ii  of 
till  pailii.id  Will'  ile'igeiili  i|  Ib-t.i  I  Ih'ta  2  tiauiiua  I*  and  ''alpha 
1  III  n-  ii  I'l  id'  iii  e  that  tin-  1' .'tipha  detiTtor  inalfunetiom-d  shurtiy  after 
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Table  I 

Iiwtrumpnl  charBcteristicB 


Iiwtni- 

ment 


[  Field  of  view 

I  M  I  AA^Q 


(deR.)  (cm*»ter) 


liCcintil-  i  Shielding; 

lAtor  I  Window ; 
(g/cm*)  I  Sides 


l>tecD>r  t.lireshold 
(MeV) 

yep 


! 

i 

3  nig/em*  \ 

A 

[ 

o.r, 

2.0 

Beta- 1 

1 

0.3 

Al;  1 

! 

90 

1.8  X 10-* : 

"B* 

>  3  g/cni"  ' 

E 

* 

i  i 

1.0 

’  8 

i 

i 

Plastic 

Al 

1 

Beta-2 

1 

C» 

w 

0.05| 

1.4 

i 

2  g/cm* 

A* 

0.5 

4.6 

40. -L 

j 

4.0 

Al; 

Ganuna 
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1.8  i 

>  4  g/cm* 

B 

* 

1.0 

5,0 

40  + 

! 

CsI(Tl) 

M 

C 

w 

0.04 

4.0 

40. 

70  rTig/em'' 

A* 

1.0 

1.0 

6.2 

12.5 

Al; 

P 

6 

1.7  X  10-* 

2.i  g.'cni* 

B 

* 

2.0 

2,0 

6.5 

CsKTl) 

Fe 

C 

4 

4.0 

4.0 

7.6 

D 

4 

8.0 

8.0 

10.5 

i 

f: 

* 

20. 

20. 

21. 

\ 

F 

* 

30 

30. 

31. 

i 

1 

G 

50.  I 

50. 

,51. 

i 

j 

H 

* 

;o. 

70. 

70. 

J 

\\ 

0.012 

0.3 

6.0 

Kty: 


4 .  Flux  nwaiiuivment. 

H' .  Kat*?  of  energy  dep'sition  nn>a8urerfient. 

A  thru  J ;  I'fuwinel  idenldioaf  :‘>n. 


launoii  Till  mont  feature.-i  of  the  n-inamitit;  tour  iuHtrutnents  are 

sumnarizpil  in  'able  i. 

It  will  In'  idisprvisl  that  deteit  irs  Ib’ta  1,  lh‘ta  2  and  (larntna  have 
extremely  wide  HehN  i>f  view  f«*i  rhargeii  p»rlie!eH  (.ieteetor  (Jainma  lias 
eaiM'iitiully  a  Iv  t'uhi  fv.f  gamma  ratliatioiO  Tlie-M-  instruments  were  iitiliyeii 
in  deDTiLsf'ng  the  gnta**  aap'-ets  nt  the  artifi'-.al  ‘adts  IVteetor  I*  •'m  an 
ertreniely  nai-Tow  held  of  view  ,  and  in  •■•sw  ntially  a  'ditss  t lona)  diflereiiti,.! 
intensitv  ”  detec  tor  This  det4*etor  w  as  u.s*mI  in  .d>tai"ii.g  |nt!  h  angle  distnbu 
tions.  Output  ehannrU  imwt  u*wsl  for  the  data  i>r»*a*’nt*’d  m  tins  iiaj.s-r  an- 
markrst  with  an  tt.'*eriak 

Al!  fleleetor*  tduamed  energy  s|w<-tr»  infi>rmation  by  rneafis  ■•}  pnls^i 
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height  diRcriiiiitiAtion  performed  in  the  variottR  electronic  logics  of  each 
output  channel,  'flireshold  diHcrimination  was  utilized,  hence  the  output 
counting  rates  an*  integral  counting  rates.  Ail  counting  rateraeters  and 
photomultiplier  current  meters  had  logarithmic  responses  in  order  to  achieve 
wide  dynamic  range.  Time  constants  of  the  circuits,  in  the  ranges  of  counting 
rate  cncount<  n  d,  wen*  less  than  I  sec. 

Discrimination  among  types  of  particles  (electrons,  protons  or  gammas) 
is  achieved  by  a  comparison  of  the  observed  counting  rates  and  the  ratM 
of  energy  dcp»-)sition  in  several  channels,  and  by  roll -modulation  observations. 
In  artificial  radiation  belts  produced  by  high  altitude  nuclear  detonations, 
the  detccterl  partic'les  are  predominantly  electnrrw. 


3.  Presentation  of  preliminary  data;  diacuaaion 

Figs.  3,  4  and  .‘i  indicate  some  of  the  gross  aspects  of  the  observed  artificial 
radiation  belts.  The  .lata  ser%'e  to  give  a  semi-quantitative  picture  of  the 
changes  occurring  in  certain  rt'gions  of  the  magnetosphere  at  that  time. 


In  tig  .4  the  utuMirT«'tt<*d  js-ak  counting  raU*  of  d'.t<*ctor  (>amma  is  given 
a.s  H  tunction  of  rln  /.  shell  for  different  tinn*s  of  olwaTvation.  The  (S'ak 
(oujitnig  rate  is  rh*  nmxiniuni  ctiumnjg  rate  ineaoure<l  during  s  tumble 
jsTusl  of  the  Iht'  III  «  »  wh<-n'  dweiling  time  m  a  particular  /.-Mhell 

suthdcntly  iojig  an  average  js>ak  value  was  ns«*d  Tile  cMunting  raU* 
w  fhu*  a  inea.HUir  o|  hut  not  iie<s'«iou'ily  linearly  ptojsirtiona!  to.  fo.'  all 
H  /,  vsliies  the  omnj  directional  Hu*.  .\ll  data  were  taken  for  an  interval 
of  magortic  field  iiit»'n,sitv  H  l>♦’t»ee;^  0  170  and  0  20.^  gauss.  The  disagree- 
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mcnt  in  counting  rates  shown  in  figs.  3  and  4  for  Zr-shcUs  lcs.s  than  1.8  arc 
probably  due  to  this  lat^e  B  interval.  The  pai  ticles  observed  by  this  detector 
are  deduced  to  lie  elcclron.s  rather  than  gamma  radiation  since:  (a)  .strong 
roll-modulation  is  apparent  in  the  raw  data;  (b)  channels  A  and  B  read 
very  much  the  .same;  and  (c)  the  rate  of  energy  depo.sition  indicatei.1  in 
channel  C  is  con.sistent  with  the  measured  flux  and  with  each  particle 
deixisiting  a  few  MeV’  of  energy  in  the  scintillator.  IVotons  of  energy  greater 
than  40  MeV  ore  ruled  out  by  the  nb.sence  of  a  signal  in  channel  F  and  by 
the  low'  counting  rates  in  channeLs  D  and  E  of  detector  P.  In  the  preparation 
of  fig.  3,  time  corrections  were  applied  to  cphcmeri.s  data  to  make  all  peaks 
coincide  at  L  1.0.  This  action  wa.s  deemed  necessary  by  the  following 
observations.  In  each  complete  satellite  orbit,  there  were  three  distinct 
regions  in  which  extremely  high  fluxes  were  encountered.  In  some  orbits, 
ephemcris  data  lerl  to  a  unique  value  of  L  as  1.9  for  each  of  these  regions. 
On  other  orbits,  three  different  L  values  were  deduet-d.  This  shift  in  L  values 
occurred  in  some  cases  even  for  consecutive  orbits.  In  all  of  the.se  latter 
ca-ses,  a  unique  time  shift  for  all  data  of  an  orbit.  u.sually  of  the  order  of 
10-30  see  (although  sometimes  in  the  hundreds  of  seconds),  led  to  the  same 
L  1.9  as.signment  for  all  three  regions.  If  this  conxetion  is  not  made, 
unrealistic  models  must  be  invoked  to  explain  llic  apparent  violent  vertical 
shifts  of  an  entire  radiation  Ixdt. 


Fig.  3  illustrates  qualitatively  the  existence  on  27  October  of  a  radiation 
belt  between  L  values  of  1.2  and  1.7,  and  another  around  1.9.  The  creation 
of  a  new  belt  around  L  1.9  on  28  October  is  evident.  There  is  some  evidence 
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of  further  injection  into  the  L  ^  1.9  belt  on  1  November.  The«e  dates  can 
be  correlated  with  known  times  of  nuclear  detonations.  All  new  artificial 
bcits  ])cak  around  £<==1.9,  and  all  exhibit  tails  extending  well  into  higher 
L  regions. 

In  fig.  4.  the  nncorrected  peak  rate  of  energy  deposition  in  detectors 
Ik'tn  1  and  2  is  given  as  a  function  of  L.  for  the  same  times  as  fig.  3.  Com- 
ni(  Ids  made  concerning  data  handling  in  the  preparation  of  fig.  3  hold  for 
th(  data  also.  The  identification  of  the  particles  as  predominantly  electrons 
rath<'r  than  gammas  is  inferred  from:  (a)  the  roll-modulation;  (b)  the 
extremely  low  efficiency  of  the  detector  for  gammas  coupled  with  the  high 
counting  rates  obaer\‘ed  in  channels  A  and  B;  and  (c)  the  consistency  of 
the  latter  counting  rates  with  the  counting  rate  of  detector  Gamma  channels 
A  and  B.  The  identification  of  the  particles  as  electrons  rather  than  protons 
is  implied  again  by  the  physics  of  nuclear  detonations,  at  least  in  the  regions 
of  the  newly  created  belts.  This  latter  statement  requires  some  qualification. 
For  L  values  less  than  about  1.7.  channels  A  and  B  of  detector  P  read 
essc'ntially  the  same.  This  indicates  the  existence  of  some  protons  in  these 
regions. 


Fig.  5.  Iswinlensity  contours  for  high  energy  elcctowis  (>  4.fi  MeV),  unoorrected  rates. 

Fig.  4  lcave.s  no  doubt  about  the  creation  of  new  artificial  radiation  belts 
on  28  October,  and  again  on  1  November,  in  the  L  9  shell.  An  addi¬ 
tional  interesting  observation  is  that  the  decay  of  the  L  fv  9  belt  between 
28  and  31  October  and  again  between  1  and  4  November  is  much  more 
pronotiuccd  than  in  fig.  3.  This  suggests  a  relative  hardening  of  the  radiation 
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in  that  belt -an  effect  not  inconsiatcMit  with  tlnnrvir'g  through  the  proccHw 
of  scattering  in  *he  atmo*^j)here  at  hi-ih  B  values. 

Fig.  5  gives  isff-intensity  ct>ntoiirs  for  electron.s  of  energy  greater  than 
4.r)  MeV  !neasun*(l  by  «lete’lf»r  (lamma.  The  iipjier  half  of  the  figure  gives 
contours  observed  on  27  ()ctofK>r:  th"  lower  half,  ‘hose  observf'tl  on  the 
281h.  The  contours  wen*  obtained  I)y  satii(ning  jK'uk  i  i.unting  lates  throughout 
the  region.  This  figure  indicates  the  “lajundaries”  for  the  artiKcial  Inlt 
created  on  the  28th  and  those  previ«»u.sly  in  existence,  accnnlited  to  an 
earlier  Soviet  .shot,  atul  to  an  AmerieHii  sfiot  of  the  summer  of  1962  l:i,  4]. 
Comment  on  fig,  5  will  l)e  curtailed  since  ‘‘'ttcr  data  pertaining  to  iso- 
inten.sity  contours  have  .since  l)eeii  obtained  and  will  be  discussed  later. 

The  remaining  data  jiresentt'd  in  what  follows  are  all  concerned  with 
pitch  angle  distributions  and  derived  jiarameters  .such  as  omni-direetioual 
flux.  The  ]iitch  angle  (Ustribution.s  were  calculated  from  Iht^  channel  A 
output  of  the  nanow  (alKUit  6°)  field  of  vit  w  dett'ctor  (detector  P)  coupled 
with  satcilitp  a.spect  relative  to  the  tnagnetic  field  vector  a.s  derived  from 
the  magnetometer  signals. 

Fig.  6  shows,  as  a  function  of  time,  a  typical  portion  of  the  raw  data 
from  detector  P  and  the  corresjronding  calculated  va'ue  of  the  cosine  ef  the 
pitch  angle  0.  If  we  acct'pt  these  pitch  angle  curves  at  face  \aliie,  we  would 
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Then*  is  little  doubt  that  the  niagnctometer  data  are  in  error.  The  angular 
di-screpancy  is  always  gnaUr  in  low  if-fiold  regions.  Fu^thermon^  if  we 
plot  the  difference  in  absoluu*  values  iK'twwn  the  magnetic  field  strength 
given  in  the  ephemeris  (calculated  from  the  Jensen  and  Cain  model)  and 
the  magnetic  field  strength  derived  from  the  magnetometers,  we  obtain  an 
oscillatory  function  with  a  time  behavior  characteristic  of  the  pitch -and -roll 
behavior  of  the  satellite.  Such  an  apparent  magnetometer  error  may  In* 
due  to  the  rc.sidual  magnetic  field  of  the  satellite  or  to  magnetometer  calibra¬ 
tion  drifts.  Insufficient  magnetomder  data  have  been  analyzed,  »o  far, 
to  enable  a  definite  coiiclasion  in  thif  rcsjiect,  although  the  latter  type  of 
error  sc'ems  to  Ik?  indicatc’d.  The  angular  disciepancy  does  not  seem  to  be 
due  to  timing  errors  since  in  the  same  ortut  the  angular  shift  required  to 
bring  about  data  con.si.stency  is  just  apt  to  be  positive  as  negative.  We 
have,  therefore,  in  our  treatment  of  the  data,  shifted  the  aspect  function 
bv  the  amount  required  to  bring  about  data  consistency.  The  magnitude 
ot  such  a  sliift  seldom  exceeds  6''  and.  a.s  stated,  tends  t<)  lie  largr  r  in  the 
low  jS-field  region-s  of  the  orbit.  The  counting  rate  curve  in  fig.  C  gives  inany 
examples  of  the  ‘Mi-callcd  “butterfly'  distribution,  which  is  characteristic 
of  radiation  b  ’.ts  generated  by  the  injection  of  particles  off  the  equatorial 
plane. 

e 
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Pm.  7.  I’lti'h  ilistriiiuiiorui;  <>1  K  -  I  M»*V,  100  ^-1200  I'T.  28  IXt.  '62; 

*’  18U  1.  2  212 

Fig  7  gives  pitch  angle  distriiiui-'ms  obtained  on  28  Octolier  in  the 
L  2.14  to  2  .12  shell  at  thre*<  different  fi  values.  The  ord  aate  is  the  raw 
.uOput  signal  vrom  ehannei  A  of  detector  1*  and  is  very  ne.irly  prii|M>rtioniiI 
to  the  logaritl  m  of  the  parti  *e  counting  rate,  Kirors  in  individual  data 
}Mimt  signals  amount  to  errors  in  counting  rate  of  the  order  of  .  10  ‘’g 
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(except  in  oases  of  telemetry  noise).  Smoothed  curves  are  probably  dependable 
to  within  5  %.  All  curves  in  the  vicinity  of  the  low  pitch  angle  cut-off  are 
somewhat  questionable  however,  not  only  because  of  the  6®  field  of  view  of 
the  detector  P,  but  moi-e  important  because  of  errors  in  the  determination  of 
the  value  of  the  cosine  of  the  pitch  angle.  The  abscissa  is  the  absolute  value  of 
the  cosine  of  the  pitch  angle,  cos  0,  plotted  from  |co3  9|  =  1  to  [cos  0|  =»  0. 
80  that  pitch  angles  read  from  0  to  90®,  left  to  right,  on  a  non-linear  scale. 
Data  points  obtained  for  90®  <  0  <  1 80“  have  been  reflected  about  90®  before 
plotting.  This  amount.s  to  assuming  fore-and-aft  symmetry  with  respect  to 
a  jB-field  line.  ThLs  assumption  appears  to  bo  valid  (at  long  enough  times 
after  a  nuclear  detonation)  in  that  the  data  points  for  9O<0<18O  plotted 
separately  (when  such  points  ore  available)  do  indeed  define  the  same  curve, 
with  fluctuations,  as  those  for  O<0<9O®. 

We  have  adopted  the  convention  of  plotting  angular  distributions  as  a 
function  of  the  cosine  for  two  reasons:  first,  solution  of  the  aspect  problem 
leads  directly  to  the  determination  of  cos  6;  second,  and  most  important, 
it  can  be  shown  that  the  area  \mder  such  curves  (when  the  ordinate  is  plotted 
linearly  as  detector  counting  rate)  is  proportional  to  the  omni-directional 
flux. 

The  curv'es  shown  obey  the  w’ell-known  transformation  law’  for  directional 
differential  intcasities  [5]: 

A  (flii  (02,  Bzy  L). 

ThLs  is  a  direct  consequence  of  Liouville’s  Theorem.  The  intensity  I  has  the 
units  of  particles/sec- cm* -ster  and  is  linearly  proportional  to  the  counting 
rate  seen  by  a  “differential”  detector  such  as  detector  P.  The  subscripts  1 
and  2  refer  to  the  angular  distributions  in  different  regions  of  the  same 
i-shell.  The  pitch  angles  0i  and  02  are  related  by 

sin*  01  sin*  02 

~Br~~Br’ 

This  results  from  the  adiabatic  invariance  of  the  magnetic  moment.  The 
angular  distributions  at  different  B  values  can  be  transformed  to  partial 
angular  distributions  in  the  equatorial  plane  by  using  these  equations. 
Fig.  8  shows  the  results  of  such  a  transformation;  the  equatorial  value  of 
B  is  denoted  by  the  subscript  0.  As  is  apparent,  all  data  points  lie  on  a  single 
curve,  “universal”  for  this  shell.  The  distribution  of  fig.  8  goes  through 
a  maximum;  the  peak  counting  rate  is  about  1000  counts/scc. 

The  authors  of  this  paper  have  adopted  as  a  standard  practice,  in  dealing 
with  pitch  angle  distributions,  tiie  immediate  transformation  of  all  distri- 
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Fi^.  8.  Einiatorial  pitcli  angle  distribution;  electrons  of  £  >  1  MeV ;  1000-1200  UT- 

28  Oct.  '62;  2.189  <  L  ^  2.212. 


butioii.s  obtaiiucl  to  the  equatorial  plane  of  the  ahell  concerned.  All  datA 
pre.sented  below  have  In^en  so  treated.  With  few  exception.^,  all  such  equa¬ 
torial  di8tribution.<»  arc  tlic  result  of  at  least  two  separate  (generally  three) 
data  acquisitions  within  the  shell,  usually  from  the  same  orbit.  Other  than 
correctioits  previously  mentioned,  no  others  were  applied  in  deriving  the 
equatorial  distril)utions.  This  procedure  leads  to  incomplete  equatorial 
distribution  functions,  since  a  distribution  (even  if  complete)  taken  at  a 
field  value  H  in  a  given  shell  can  only  define  the  equatorial  distribution  up 
to  a  “cut-off”  angle  given  by 

sin*  =  BcIB. 

That  Is,  no  information  can  be  obtained  on  the  pitch  angle  uistributions 
of  those  particles  which  miiror  within  the  shell  between  the  region  of  B 
(or  its  conjugate  point)  and  the  equatorial  plane  at  Bo.  A  few  distnbutions 
arc  incomi)lctc  because  of  incomplete  ])itch  angle  coverage  by  the  satellite 
pitch  and  roll.  A  com])letc  distribution  at  some  point  Bi  will  always  lead 
to  a  complete  distribution  at  all  jwints  in  the  same  shell  w’hich  have  a 
B>Bi,  of  course. 

We  have  found  the  aforementioned  procedure  very  irseful.  Besides  limiting 
the  number  of  variables  which  must  be  considered  (to  fj,  I,  cos  6o.  L  in 
equatorial  plane)  it  facilitates  identification  of  artificial  belts,  and  allows 
checking  of  data  consistency  in  a  standard  way.  (All  di.stributions  obtained 
from  a  given  1,-shell  at  essentially  the  same  time  must  transform  to  the 
same  equatorial-plane  distribution.)  It  further  provides  a  convenient  function 
for  generating  and  determining  omni-directional  fluxes  in  regions  of  space 


1  T7 


I\'.  1:.'!  mk\.si  ki;mknts  on  tkv'  I'kii  ivvktk  i  ks  (i:»7 

iuaccfHsible  to  a  .>ai  H.U'.  'Fhis  is  acconijilislusl  liv  first  tmiisforniiiif;  the 
<listril)itti()n  from  {Hq.  L)  to  (B.  L).  and  then  intf-giatiii"  over  angle. 
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27  Oct.  <>2 


Fig.  9  gives  Oijnatorial  pitch  angle  jlistrihuthuis  for  L  shell.s  in  the  vicinity 
of  1.8,  obtained  on  27  OettilK'r.  The  d«>t4'ctor  P  calibration  has  been  in¬ 
corporated.  so  that  the  ordinatt*  is  in  units  of  actual  detector  counting  rate. 
Thi.s  can  be  converted  to  diri'ctional  dift'ercntial  intensity  (in  units  of  jiani- 
clcs'sco  cm^  stcr)  by  nniltiplieatiop  by  5  9  <  102.  \|)  curve.s  are  smoothed, 
but  the  scatter  of  data  points  is  about  the  same  a.-!  .shown  in  ng.  H.  All  the 
angular  distributions  exhibit  what  may  In-  cailcd  a  ‘‘jiredetonation’’  shajie. 
This  distribution  favors  pitch  angles  clo.ser  to  90’  and  iir's  a  low  angle  cut-off 
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Fig.  10  also  shows  diatributiojis  for  L-shells  in  the  vicinity  of  l.S.  but 
obtained  a  day  later  on  2S  October.  A  comparison  with  fig.  9  indicates  t)ic 
injection  of  particle.s  into  the  shells  in  this  vicinity  which  occurretl  some¬ 
time  Ix-tween  the  two  data  acquisitions.  The  major  injection  is  localized 
within  the  to  1.835  shell,  although  some  injection  into  both  the 

L  as  1.80  and  the  L  sn  2.20  shell  is  also  obvious.  (The  curve  for  L  2.20 
is  the  same  as  that  shown  in  fig.  8.)  The  time  can  be  correlated  with  the 
detonation  of  a  Soviet  nuclear  device.  All  distributions  ilhwlrate  the  “butter¬ 
fly”  sha])e  mentioned  earlier,  being  most  pronounced  for  L-^  l.EJS  to  1.835. 

Simiile  qualitative  arguments  can  be  usetl  to  explain  the  “butterfly” 
shajic,  and  it  may  be  well  to  jxiint  them  out  here.  Assume  point  injection 
of  electrons  following  the  detonation  of  a  nuclear  device  at  an  extremely 
high  altitude  well  above  the  tangible  atmosphere.  Assume  that  the  atmos]ihere 
can  lie  represented  by  a  slab  and  that  any  particle  striking  the  top  of  this 
slab  is  lost.  Let  Bi,  Bi  and  Ba,  be  the  magnetic  field  strength-s  at  the  top 
of  the  slab,  at  the  point  of  injection  and  in  the  equatorial  plane,  respectively, 
of  the  shell  of  injection.  Then  the  equatorial  pitch  angle  distribution  of 
electrons  will  lie  characterized  by  a  iow  angle  cut-off,  6t,  given  by 


sin*  &r 

(i.e.,  no  particle  mirrors  below  the  top  of  the  slab)  and  a  high  angle  out-off, 
01,  given  by 

sin*  01  Bo/Bi 


(i.e.,  no  particle  mirrors  between  Bi  or  its  conjugate  point  and  B©)-  The 
sharpness  of  the  low  angle  cut-off  is  smeared  out  since  the  atmosphere  is 
not  an  ideal  slab  with  a  definite  lop.  The  sharpness  of  the  high  angle  cut-off 
is  smeared  out  since  the  injection  does  not  occur  at  a  point  but  has  very 
significant  extent  in  space  due  to  the  expansion  and  possible  trapping  of 
the  fission  debris.  Radiation  previously  trapped  in  the  shell  with  a  more 
“normal”  distribution,  (pitch  angles  nearer  90”  favored)  completes  the 
smearing  out  of  the  high  angle  cut-off  to  90®. 

Consider  next  a  detonation  at  a  low  altitude  within  the  atmosphere.  The 
only  major  contribution  to  the  radiation  populating  the  belt  for  any 
appreciable  iieriod  of  time  will  be  that  due  to  the  ^-decay  of  fission  debris 
isotopes  which  have  been  transported  out  of  the  dense  region  of  the  atmos¬ 
phere.  In  this  case  we  obtain  a  distribution  of  the  same  "butterfly”  shape 
as  before. 

It  will  be  noted  in  fig.  10  that  the  distribution  for  i/»  1.818  to  1.836 
peaks  at  about  37®;  another  distribution  obtained  for  this  shell  earlier  on 
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28  October  j»cRk.s  around  32®.  I’articlefl  of  pitch  angle  S.'S®  mirror  at  a  magnetic 
latitude  of  30°  and  mean  altitude  of  about  2300  km.  The  direct  pliysicnl 
sigtuficance  of  this  is  not  obvioua  since  the  mea-surementa  wen?  |K’rformod 
at  relatively  late  times  after  the  <letoj\ation.  The  measiiremcnt  of  the  angtilar 
distribution  of  tlie  injected  particles  at  early  times  after  detonation  may 
prc.sumably  yield  information  on  the  ittjection  ntechanism. 


Fig.  11.  Oinni-diroctional  flux;  electrons  of  £  >  1  MoV’. 


Fig.  11  is  derived  from  the  data  presented  in  figs.  9,  10  and  12,  and  gives 
the  omni-dircctionol  flux  as  a  function  of  X-shcll  before  the  dctonatioit  for 
BjBo  equal  to  1.58,  and  after  the  detonation  for  BjBo  equal  to  1.58  and 
2.35.  BjBo  equal  to  unity  cannot  l)c  used  .since  the  equatorial  pitch  angle 
distribution-s  are  incomplete,  and  it  is  desirable  to  avoid  extrapolation.  The 
choice  of  a  cut  through  L-shells  at  a  constant  BjBo  is  made  bccau.se  for 
identical  equatorial  ])itch  angle  distribution.s  in  various  shells,  sucli  a  cut 
yiekljs  also  identical  pitch  angle  distributions.  This  choice  thus  seems  the 
natural  basis  for  comparison  when  it  is  desired  to  show  the  effects  of  shell 
numl»cr  only.  Fig.  H  is  the  “coiTcct”  version  of  figs.  3  and  4  showm  earlier. 
Again,  tiie  belt  created  on  28  October  is  clearly  evident,  v’ith  a  peak  omni¬ 
directional  flux  in  Z/=1.83. 

Fig.  1 1  shows  a  feature  not  evident  in  figs.  3  or  4  although  suggested  by 
the  structure  beyond  L  —  2.0.  Differences  between  these  representations 
should  not  be  unexpected,  however,  since  the  figim»  rcjirescnt  different 
monies  of  cutting  through  L-shclls  and  the  figures  refer  to  different  threshold 
energies.  Moreover,  the  narrow  field  detector  yields  a  more  accurate  measure¬ 
ment  of  omni-directional  flux.  The  1.58  curve  .shows  a  .slot  at  1.05 

and  a  subsequent  rise  in  counting  rate  which  presumably  goes  through  a 
j>cak  beyond  L  —  2.05.  The  fi/J5o'=2.35  curve  shows  clearly  the  existence  of 


HO 


L.  KATZ  ET  AL. 


(tOl) 


II V.  12 


a  secoad  belt,  peaked  at  L  2.10,  aiid  somewhat  broader  than  that  at 
Z(=1.83.  Both  belts  were  apparently  generated  by  the  same  nuclear 
detonation,  judging  from  the  1.58  ctiree  taken  on  the  27th  of  October. 

Fig.  12  givc.s  iso-intensity  contours,  on  a  dipole  field,  for  28  October. 
The  number  defining  the  contours  arc  omni-directional  fluxes  derived  by 
first  transforming  the  equatorial  pitch  angle  distribution  of  each  L-shcll  to 
various  B  values  in  that  shell,  and  then  integrating  over  pitch  angle  at  each 
B  value.  Elach  ntimber  used  is  probably  good  to  within  10  %  considering 
only  .statistical  accuracy.  Larger  systematic  errors  previously  discussed  may 
exist  in  the  determination  of  position.  The  omni-directional  flux  J  at  a 
j>oint  {B,  L)  is  given  explicitly  by 

0 

J(2?,  X)  =  4jt  /  /(cos  0,  B,  L)  d{coa  fl). 


Fig.  12  is  the  “correct"  version  of  fig.  5.  One  feature  that  at  first  glance 
may  seem  surprising  —  namely  a  region  of  maximum  intensity  in  the  1.83 
shell  off  the  equator,  around  .6  =  0.08  — is  a  natural  consequence  of  the 
“butterfly”  shajx!  of  the  distribution.  In  fact,  any  distribution  that  presents 
a  ^-cak  at  an  angle  different  from  90®  in  the  equatorial  plane  will  lend  to  a 
maximum  omni-directional  flux  at  some  high  B  value  in  the  shell.  The 


Fig.  12.  I»o-intrn«ity  coiitonrs;  electrons  of  S  >  1  MoV;  1000-1200  UT.  28  Oct.  ’62. 

contours  in  the  region  defined  by  1.2<L<1.6,  and  0.26<B<0.16  are 
shown  as  “uncertain”  since  there  are  indications  from  our  data  that  these 
arc  due  to  an  adntixture  of  both  protons  and  electrons.  Insufficient  analysis 
has  l)eon  performed  for  an  accurate  determination  of  the  relative  abundance 
of  each  particle.  These  arc  probably  contours  of  the  natural  inner  belt. 
Omni-directional  fluxc.s  observed  arc  not  inconsistent  with  an  proton 
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integral  spectrum ,  with  E  =  ^  MeV,  where  fluxes  observed  from  prot^uis  (»f 
E>40  MeV  by  Van  Allen  [6]  arc  used  to  determine  the  scaling  factor 
(about  6  <  lot). 

Fig.  12  indicates  dearly  the  existence  of  two  well  defined  belts  on  2h 
Oc  ober.  A  Soviet  detonation  was  known  to  occur  in  the  vicinity  of  L  1.H3 
on  this  date.  The  In-lt  at  L  1.83  i.s  the  shsrp«‘r  of  the  two,  and  the  onr  most 
easily  recognixable.  The  shell  at  L  2.1  i.s  neverthele.s.s,  also  qjiite  is'al.  and 
new,  as  data  from  27  October  (see  e.g.  figs.  7  and  ft)  show  low  particl- 
counting  rates  in  this  shell. 


•» 


13.  F>)U<ilori«l  f)!teh  an(;ie  «{iHfr'*)ui;'>ns;  L  1,83. 

Figs.  i3  through  15  |N>rtain  to  the  time  In-havior  of  the  L  1.83  radiation 
IsOt  iK'tween  28  and  30  Oetol>(*r.  The  ordinates  of  figs.  13  and  14  are 
dcfiendable  to  within  5  except.  ii<  usual,  in  the  vicinity  of  tlie  cutotl; 
the  ordinate  of  fig.  1,'  is  probably  good  to  7  Fig.  13  givc.s  the  c(pritorml 
distribution  at  L  1.83  at  thre«'  different  tinie.s  after  the  event,  and  on  i  he 
day  Indore  the  event  The  most  obvious  features  aix'  ttie  sfnft  of  the  jk  uk 
trward.s  !I0  and  the  decrease  in  the  js-ak  (a  tendeney  toward.s  isotropy)  ni 
ijualitative  agn'cment  with  the  more  pronoutusnl  etb  ct  e.f  atmospheric 
Rcattering  on  parti'  le!-  of  lower  pitih  angle  The  cpiestiou  of  whethei  these 
di'cay  rat«“s  are  consi.steu?  with  removal  t>v  scattering  on  .atniosphei ic  particles 
alone  is  Inniig  stmlied  and  is  obviously  of  geophysical  signiii.  ance 

Fig.  14  IS  a  plot  of  inti  nsities  in  the  /.  1  83  shell  versus  time  for  diffenmt 

pitch  angles  ('hang«‘s  an-  s<-cn  to  occur  most  tapidly  at  earlv  times  partn  ie 
inten.sit ICS  at  the  smaller  pitch  angle-  temj  to  <-hangc  tlic  faste.st  .  and 
distributions  at  lat«'  times  tend  to  apj'roaeh  isotropy  (a.s  indicated  bs  the 
groupings  of  jMiints  at  and  5.n  hr)  riic  rxTi  melv  fast  dM'a>  of  the  intensity 
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r»t  pilch  angle  32".  which  mirrors  at  a  mean  altitude  of  190')  km,  is  not  too 
.surjjrising.  The  initial  fast  decay  of  the  intensity  at  pitch  angle  43“,  which 
mirrors  at  a  mean  altitude  of  3100  km.  i.s  suriirising.  If  the  data  are  believed, 
this  suggests  a  reinftvui  mechanism  other  than  atmospheric  scattering  for 
this  L-shell  ngic'n  although  a  definite  conclusion  mu.st  Ik-  reserved  until 
atnuispheric  scattering  calculations  are  {K-rformed.  It  would  also  lie  well 
to  look  closT'ly  at  pitch  angle  distributions  of  radiation  in  unj>effurl)e<i  shells. 
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h'Uf  l.'i  ,  Kh-iic  Ilf  /  1  S.'t  Ix  it  iitiiiii  iliOTti.iiinI  flux,  /<  Ho  i 

and  study  th  -m  fur  li  ng  {m  i-kmIs  of  time.  Th*'  ndative  time  nule)s-nd'‘Ute 
uf  the  inti'osi  ics  at  pitch  angles  rt2  and  fOt  .  (and  higher,  as  tndicat<'<l  in 
fig.  13)  is  uiid  rstandaiiii'.  There  n  ay  *’ven  Ik*  a  tendency  for  int<‘nsiti»‘s  at 
verv  iargi'  pitch  angles  tu  incn-aKi-  with  tune  (as  suggeste<1  bv  the  i»i'  curve: 
ati  effi'Ct  i-uni  isteiit  with  tin-  llieun-fical  caleulatiuiis  uf  Walt  IT]  showing 
upwanl  UK  well  as  duwtiward  diffusion  of  mirror  {sur.ts  of  partich's  scrttt«rrd 
in  till'  atnioH  plierc 
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Fig.  15  is  a  replot  of  the  data  of  fig.  13,  giving  omni-directional  llu.Te.s 
in  the  L=‘  1.83  l»elt  as  a  function  of  time.  The  decay  of  the  radiation  com¬ 
ponent  injected  into  this  shell  on  28  October  i.s  not  inconsistent,  for  the 
three  day  interval  of  observation,  with  an  e.xponential  decay  half-life  of 
approximately  75  hr.  However,  the  decay  rate  is  certainly  very  deiKuidcnt 
on  pitch  angle  distribution.  Since  all  distributions  tend  ultimately  towanl 
the  “normal”  shape,  with  pitch  angles  towards  90“  favored,  the  decay  rate 
should  certainly  decrease  with  time.  ComparLson  of  figs.  14  and  15  shows 
that  the  drastic  effects  indicated  in  fig.  14  tend  to  be  masked  by  the  integra¬ 
tion  over  pitch  angle,  and  strongly  implie.s  that  curves  like  fig.  14,  derived 
from  pitch  angle  distributions  are  most  iLseful  in  time  behavior  studies. 

4.  Closing  remarks 

The  possibility  of  the  creation  of  double  belts  by  a  single  nuclear  detona¬ 
tion  is  indicated  by  the  data.  We  have  as  an  example,  in  the  iircsent  ca.se. 
the  generation  of  two  well  defined  belts,  at  L  1.83  and  Z.  ^  2.1,  by  the 
burst  of  28  October. 

Mechanisms  for  electron  removal  other  than  atmospheric  scattering  may 
be  at  work.  In  any  event,  the  most  jirofitable  means  of  studying  the  time 
behavior  of  artificial  belts  seems  to  be  by  futeasity  at  constant  pitch  angle 
versus  time  data. 

Information  on  omni-directional  flux  and  iso-intensity  contour  maps 
alone  have  limited  usefulness  as  a  means  of  studying  the  dynamics  of  radia¬ 
tion  belts.  Data  of  prime  importance  are  pitch  angle  distributions  obtained 
at  various  L  values  in  the  equatorial  plane.  These  can  be  used  to  obtain 
the  omni-directional  flux  and  iso-intensity  contoiurs  as  well. 

Thb  rejjort  may  serv’c  to  illustrate  the  methods  of  da^  handling  used 
by  the  authors,  and  may  give  some  indication  of  the  nature  of  the  information 
to  be  learned  from  satellite  data  analysis,  in  particular  from  pitch  angle 
distributions.  Work  will  continue  in  all  asjiects  of  data  reduction.  Of  particular 
interest  will  be  data  obtained  subsequent  to  the  latest  observed  nuclear 
detonation  which  wall  allow  a  long  time  study  of  the  artificial  belt  created 
by  that  detonation  without  complication  by  further  injection.  Data  ]X!r- 
taining  to  energy  sjxictra,  so  far  almost  completely  untouched,  will  l>o 
analyzed.  Derivation  of  iso-intensity  contours,  using  the  methods  of  pitch 
angle  transformation  described  in  this  report,  are  in  progress  for  various 
times  of  interest.  Attempts  will  bo  made  to  calculate  the  effects  of  atmospheric 
scattering  on  pitch  angle  distribution.s,  and  to  check  the  consistency  of 
the  results  of  such  calculations  with  the  experimentally  observed  data. 
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A  minor  program  to  detf'ct  systeroatio  errors  in  the  magnetometer  calibrations 
will  also  be  attempted  in  hojws  of  improving  aspect  determination. 
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